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SYMBOLS 

Exchanger  total  heat  transfer  area  on  one  side,  ft2 
Exchanger  minimum  free-flow  area,  ft2 
Exchanger  total  frontal  area,  ft2 
Plate  thickness,  in 
Plate  spacing,  in 

Flow  stream  capacity  rate,  BTU/(hr  °F) 

Jet  contraction-area  ratio,  dimensionless 
Capacity  rate  ratio 

Specific  heat  at  constant  pressure,  BTU/(lbm  °F) 

Air-side  header  diameter,  ft 
Mean  friction  factor,  (Equation  32) 

Fuel -air  ratio 

Exchanger  flow-stream  mass  velocity,  1 bm/ (hr  ft2) 

Proportionality  factor  in  Newton's  second  law,  lbm  ft/(lbf  sec2) 

Unit  conductance  for  thermal -convection  heat  transfer, 

BTU/(hr  ft2  °F) 

Colburn  factor  *  NSj  NpR 
Bend  loss  coefficient 

Contraction  loss  coefficient  for  flow  at  heat  exchanger  entrance, 
dimensionless 

Momentum  flux  correction  factor,  dimensionless 

Expansion  loss  coefficient  for  flow  at  heat  exchanger  exit,  dimensionless 

Constant  used  In  Equation  14 

Constant  used  In  Equation  15 

Unit  thermal  conductivity,  BTU/(hr  ft2  °F/ft) 
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L  -  Heat  exchanger  counter-flow  length,  ft 
L  -  Gas-side  header  length,  ft 

Ln  -  Heat  exchanger  non-flow  length,  ft 
i  -  Fin  length  from  root  to  center,  ft 

M  -  Molecular  weight 

m  -  A  fin  effectiveness  parameter  {2h/k5)1^ 

NpR  -  Prandlt  number,  dimensionless 

NR  -  Reynolds  number,  dimensionless 

N$t  -  Stanton  number,  dimensionless 

NTU  -  Number  of  transfer  units 

P  -  Power 

Pf  -  Fractional  Power 

p  -  Pressure,  lbf /in2 

o 

q  -  Dynamic  velocity,  lbf/ft* 

Ru  -  Universal  gas  constant,  ft  lbf/(lb  mol)(°R) 
rh  -  Hydraulic  radius  *  ACL/A 

T  -  Temperature,  °F  or  °R 

U  -  Unit  overall  thermal  conductance,  BTU/(hr  ft2  °F) 

V  -  Volume,  ft3 

V  -  Velocity,  ft/sec 

v  -  Specific  volume,  ft3/lbm 

W.  -  Weight  of  air-side  fins  and  plates,  lbf 

d 

WENCL  *  Het9ht  of  enclosure,  lbf 

Wg  -  Weight  of  gas-side  fins  and  plates,  lbf 

WHD  "  Weight  of  headers,  lbf 

Wj  -  Total  weight  of  heat  exchanger,  lbf 
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w  -  Hass  flow  rate,  lbm/sec 

a  -  Ratio  of  total  transfer  area  on  one  side  of  the  exchanger  to  total  volume 
of  the  exchanger,  ft2/ft3 

0  -  Ratio  of  total  heat  transfer  area  on  one  side  of  a  plate-fin  heat 

exchanger  to  the  volume  between  the  plates  on  that  side,  ft2/ft3 
a  -  Denotes  a  difference  between  values  of  the  same  parameter 
6  -  Fin  thickness.  In 

e  -  Heat  exchanger  effectiveness,  dimensionless 

Of  -  Fin  effectiveness,  dimensionless 

n0  -  Total  surface  effectiveness,  dimensionless 
a  -  Ratio  of  free  flow  area  to  frontal  area 

u  -  Viscosity  coefficient,  lbm/(hr  ft) 

if  -  Constant  =  3.1416 

p  -  Density,  lbm/ft3 

pc  -  Density  of  core,  lbm/ft3 

%  -  Density  of  material,  lbm/ft3 

p„  -  Density  of  material  per  unit  area,  lbm/ft2 

-  Fin  weight  factor 
i|u  -  Plate  weight  factor 
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SUBSCRIPT 


a  -  air-side 

av  -  average 

b  -  bend 

c  -  core 

g  -  gas-side 

1  -  inlet 

m  -  mean  value 

max  -  maximum  value 

min  -  minimum  value 

o  -  outlet 

t  -  total 

1  -  entrance  condition 

2  -  exit  condition 

3  -  compressor  exit  condition 

4  -  turbine  Inlet  condition 

5  -  turbine  outlet  condition 

6  -  recuperator  outlet  condition 
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LIST  OF  ABBREVIATIONS 


BTU 

- 

British  thermal  unit 

COC 

- 

Control  Data  Corporation 

°F 

- 

Degrees  Fahrenheit 

ft 

- 

Feet 

hr 

- 

Hour 

In 

- 

Inch 

lbf 

- 

Pounds  force 

- 

Pounds  mass 

NTU 

- 

Number  of  Transfer  Units 

psla 
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Pound  force  per  square  inch  absolute 

°R 
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Degrees  Rankine 

sec 

_ 

Second 
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ABSTRACT 


This  report  presents  a  computer  program  for  preliminary 
design  analysis  of  counter-flow,  compact,  plate-fin  heat 
exchangers.  The  program  method  is  based  on  the  effectiveness-NTU 
relationship  analysis.  The  heat  exchanger  design  begins  with 
assumptions  for  counter-flow  length,  total  frontal  flow  area  and 
core  matrix  fin  geometry.  Using  these  constraints,  the  program 
proceeds  to  calculate  the  resulting  effectiveness  and  pressure 
drop  based  on  specified  air-side  and  gas-side  conditions.  Input 
design  requirements  include  selected  air-side  and  gas-side  fin 
types;  the  pressures,  temperatures,  and  mass  flows  of  the  air  and 
gas  streams;  fuel -air  ratio;  as  well  as,  the  maximum  air-side 
inlet  header  velocity.  Heat  exchanger  designs  may  be  generated 
based  on  four  different  fin  types  ( i . e . ,  plain,  louvered, 
strip/offset  or  wavy  fins)  over  a  varied  number  of  core 
dimensions.  Program  output  includes  inlet  and  exit  conditions  on 
air  and  gas  sides,  effectiveness,  fin  characteristics,  core  length 
and  volume,  total  frontal  flow  area,  pressure  drops,  overall 
enclosure  height,  number  of  transfer  units,  overall  weight,  and 
air-side  header  diameters  and  velocities.  This  report  presents 
the  analysis  method,  description  of  input  and  output  with  sample 
cases,  and  a  program  listing. 

ADMINISTRATIVE  INFORMATION 

This  documentation  was  accomplished  under  Work  Unit  1-2720-150,  Element 
62543N,  Task  Area  S  0340  SL039,  Task  23556.  The  project  is  a  part  of  the 
Propulsion  Technology  portion  of  the  Ships,  Submarines  and  Boats  B  ock  Program. 
The  work  was  done  in  the  Engines  Branch  of  the  Power  Systems  Division,  Propulsion 
and  Auxiliary  Systems  Department,  David  Taylor  Naval  Ship  Research  and  Development 
Center.  Program  manager  at  the  Naval  Sea  Systems  Command  was  Mr.  D.  A.  Groghan 
(SEA  05R13) . 
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INTRODUCTION 


Preliminary  design  analysis  of  compact  heat  exchangers  for  gas  turbine  engine 
applications  Involves  repetitious  calculations  varying  design  and  performance  over 
a  range  of  conditions.  In  order  to  determine  "an  optimum"  design,  thousands  of 
these  calculations  nust  be  completed.  Design  accuracy  and  details  which  are 
necessary  In  manufacturing  a  heat  exchanger  are  not  necessary  for  this  type  of 
preliminary  design  analysis.  General  and  approximate  procedures  are  sufficient  to 
yield  the  desired  design  and  performance  characteristics. 

This  report  presents  a  computer  program  for  preliminary  design  analysis  of 
counterflow,  compact,  plate-fin  heat  exchangers.  The  computer  program  Is  based  on 
the  effectiveness-NTU  relationship,  [1].  The  heat  exchanger  design  begins  with 
assumptions  for  counter-flow  length,  total  frontal  flow  area  and  core  matrix  fin 
geometry.  Using  these  contraints,  the  program  proceeds  to  calculate  the  resulting 
effectiveness  and  pressure  drop  based  on  specified  air-side  and  gas-side  Inlet 
conditions.  Input  requirements  are  air-side  and  gas-side  fin  types;  the 
pressures,  temperatures,  and  mass  flows,  of  the  air  and  gas  streams;  fuel -air 
ratio;  as  well  as,  the  air-side  inlet  header  velocity.  The  program  gives  the 
designer  the  capability  to  select  from  four  different  fin-types  (i.e.,  plain, 
louvered,  strip/offset,  and  wavy)  with  a  variety  of  surfaces  (I.e.,  fins/in  and 
fin  heights)  for  each  fin  type.  All  the  necessary  data  and  characteristics  for 
the  fin  types  are  located  in  reference  1.  The  Input  data  Includes  assumed  values 
of  core  counter-flow  length  and  frontal  area  to  be  used  In  the  computations 
described  below.  Output  consists  of  all  input  requirements,  as  well  as, 
calculated  parameters  Including  core  volume,  pressure  drops,  overall  enclosure 
height,  air  and  gas  exit  temperatures  and  pressures,  number  of  transfer  units, 
heat  exchanger  effectiveness,  overall  weight,  and  air-side  header  exit  diameter 
and  velocity. 

A  complete  description  of  Input  and  output  variables,  a  FORTRAN  IV  program 
listing,  and  discussion  of  the  analysis  method  are  Included  In  the  report.  Also, 
Included  are  examples  for  using  the  program,  and  illustrations  of  output. 
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PCTHOO  OF  ANALYSIS 


In  sizing  heat  exchangers  there  are  two  parameters  which  affect  size  and 
shape.  These  parameters  are  effectiveness  e,  and  pressure  drop  dp. 


The  objective  of  using  a  heat  exchanger  as  a  regenerator  In  gas  turbines  Is 
to  raise  the  compressor  exit  temperature  using  the  waste  heat  from  exhaust  gases; 
therefore.  Increasing  thermal  efficiency.  Raising  the  air  temperature  In  a 
regenerator  to  that  of  the  entering  gas  temperature  would  constitute  a  perfect 
heat  exchanger.  How  close  the  air-side  exit  temperature  reaches  the  entering  gas- 
side  temperature  defines  effectiveness.  In  order  to  define  heat  exchanger 
effectiveness,  a  capacity  rate  Is  used.  The  capacity  rate  Is  the  mass  flow  rate 
times  the  heat  capacity.  For  regeneration  In  gas  turbine  engine  applications,  the 
capacity  rate  Is  slightly  lower  on  the  air-side  then  on  the  gas-side  due  to 
compressor  bleeds  lowering  air-side  mass  flow  rate,  combustion  products  Increasing 
gas-side  mass  flow  rate,  and  Increasing  heat  capacity  as  temperature  rises. 
Effectiveness  Is  defined  as  [1], 


Alternatl vely. 


T.  *  air-side  exit  temperature 
®o 

Ta^  =  air-side  entering  temperature 

T  a  gas -side  exit  temperature 

"o 

Tg^  *  gas-side  entering  temperature 

C  *  capacity  rate. 


(la) 


(lb) 


Heat  transfer  can  be  related  to  a  nondlmensional  parameter,  known  as  the 
number  of  transfer  units  (NTU),  In  terms  of  [1], 


NTU  * 


A  U 


av 


'min 


(2) 
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where, 

A  *  heat  transfer  area,  ft2 

Uay  =  average  overall  heat  transfer  coefficient,  BTU/(hr  ft2  °F) 
cm1n  s  capacity  rate,  BTU/(hr  °F) 

For  a  counterflow  heat  exchanger,  which  this  computer  program  Is  based  upon,  the 
relationship  of  effectiveness  to  NTU  Is  [1], 


-  NTU  (1-  CB(n/CMx) 

£■^-4 - .  (3) 

1  -  r~  ,  -  “Tu  t‘-  WCTOx> 

max 

Equations  for  other  configurations  (I.e.,  crossflow,  parallel  flow,  etc.)  can  be 
found  in  reference  1. 

Pressure  drop  (Ap)  In  a  heat  exchanger  can  adversely  affect  the  performance 
of  the  heat  exchanger  and  the  specific  power  of  the  gas  turbine.  Pressure  drop 
can  be  made  nondlmensional  by  dividing  the  change  In  pressure  by  the  absolute 
pressure  before  the  pressure  drop  occurs.  This  nondlmensional  pressure  drop  can 
be  related  to  a  fractional  pressure  drop,  which  is  directly  related  to  a  power 
loss  In  the  gas  turbine,  regardless  of  whether  the  4P  occurs  on  the  air-side  or 
gas-side,  [2].  The  above  mentioned  statement  may  be  made  clearer  with  the  aid  of 
figure  1  and  the  following  derivations:  At  a  given  turbine  Inlet  temperature, 
power  Is  proportional  to  the  pressure  ratio  across  the  turbine  (P4/P5).  Assume 
that  an  overall  pressure  ratio  (p3/P6)  can  be  related  to  P4/P5  by  considering  the 
pressure  losses  across  the  heat  exchanger  and  burner  (dpa  *  air-side  drop,  tyg  = 
gas-side  drop,  and  ^  =  burner  drop);  then  the  recuperated  turbine  pressure  ratio 
can  be  written  In  terms  of  the  overall  pressure  ratio. 


P4  P3  -  4>b  -  4>a  P3  k  P3  P3 

P5“P^  (l+i) 


(4) 
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The  fractional  power  loss  due  to  recuperation  can  be  expressed  In  the  following 
manner  for  small  pressure  losses. 


where. 


(6a) 


(since  Ap 

a 


(6b) 


RECUPERATED 


B  -  BURNER 
C  -  COMPRESSOR 
R  -  RECUPERATOR 
T  -  TURBINE 

Q  STATION  NO. 


Figure  1  -  Simple  {Unrecuperated)  and  Recuperated  Gas  Turbine  Cycles 
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Substituting  In  Pnew  and  P^,  remembering  that  the  pressure  drops  are  small 
relative  to  the  absolute  pressures,  and  Ignoring  second  order  terms,  equation  (5) 
becomes. 


P3  Pe 


(7) 


Therefore,  the  fractional  power  loss  due  to  regeneration  Is  related  to  the 
fractional  pressure  losses  across  the  recuperator  by  equation  (7).  Since  power 
losses  In  the  cycle  can  result  from  pressure  drops  in  the  heat  exchanger.  It  is 
desirable  to  keep  pressure  losses  to  a  minimum.  This  can  be  accomplished  by 
choosing  fin  types  that  do  not  create  high  pressure  losses.  For  example,  plain 
fins  Instead  of  strip  fins,  short  flow  lengths,  or  large  frontal  areas,  which  may 
create  other  problems  that  are  not  dealt  with  In  this  manual. 


As  stated  earlier,  the  heat  exchanger  design  begins  with  assumptions  for 
counter-flow  length,  total  frontal  flow  area  and  core  matrix  fin  geometry.  Using 
these  constraints,  the  program  proceeds  to  calculate  the  resulting  effectiveness 
and  pressure  drop  based  on  specified  air-side  and  gas-side  Inlet  conditions.  The 
basic  procedure  Is  summarized  below: 

(1)  GIVEN:  Mass  flow  rates,  Inlet  pressure,  and  Inlet  temperatures  on 

both  air  and  gas  sides. 

(2)  SELECT:  Surface  characteristics  (1 . e. ,  f Ins/In  and  fin  height)  both 

sides 

(3)  CALCULATE:  Heat  transfer  and  free  flow  areas  on  both  sides. 

Equations  (9)  and  (II) 

(4)  DETERMINE:  Fluid  properties  on  both  sides 

(5)  CALCULATE:  Reynolds  number  of  both  sides,  Equation  (13) 

(6)  DETERMINE:  Stanton  number,  Colburn  factor  and  friction  factor 

(7)  CALCULATE:  Heat  transfer  coefficient  on  both  sides.  Equation  (17) 

(8)  CALCULATE:  Fin  effectiveness  on  both  sides.  Equation  (18) 

(9)  CALCULATE:  Surface  effectiveness  on  both  sides.  Equation  (21) 
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(10)  CALCULATE:  Overall  coefficient  of  heat  transfer  based  on  the  air-side 

area.  Equation  (22) 

(11)  CALCULATE:  NTU  and  exchanger  effectiveness,  Equations  (2)  and  (25) 

(12)  CALCULATE:  Inlet  and  exit  loss  coefficients,  Equations  (34)  and  (35); 

and  core  pressure  drop  on  both  sides.  Equations  (30),  (36), 
(37),  and  (38) 

This  procedure  requires  the  following  Inputs:  air  flow  rate  wa,  fuel -air  ratio 
far,  air-side  fin  surface,  air-side  entering  pressure  pla,  air-side  entering 
temperature  T^a,  gas-side  fin  surface,  gas-side  entering  pressure  p^g,  gas-side 
entering  temperature  Tlg,  a  plate  thickness  a,  and  a  thermal  conductivity  k  for 
the  fin  material  and  parting  plates.  By  specifying  the  fin  surface, 
characteristics  accompany  the  fins,  such  as,  plate  spacing  b,  hydraulic  radius 
rh,  fin  thickness  6,  ratio  of  transfer  area  to  volume  between  plates  B,  and  ratio 
of  fin  or  extended  area  to  total  area  are  determined.  The  frontal  area  Afr,  and 
the  heat  exchanger  core  volume  ¥  must  be  specified  for  the  calculations. 

HEAT  TRANSFER  AND  FREE -FLOW  AREA 

The  ratio  of  total  transfer  area  on  one  side  (i.e.,  air-side,  gas-side)  of  the 
heat  exchanger  to  the  total  volume  of  the  exchanger  is  given  by, 

“a  =  777 T  =  b  A  3+  2a  *  for  the  air'side’  <8a> 

total  a  g 

and 

A  b  B 

°g  =  7777  =  b  "  A  +  2a  *  for  the  9as_si  de*  (8b) 

Rearrangement  of  Equation  (8a)  or  (8b)  gives  the  total  heat  transfer  area  on  one 
side. 


A  3  “  *total  • 

The  ratio  of  free-flow  area  to  frontal  area  Is  defined  as. 


(9) 
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(10) 


Rearrangement  of  Equation  (10)  give  the  free-flow  area  on  one  side  as 

Ac  =  a  Afr  *  {11) 

FLUID  PROPERTIES 

An  Initial  value  for  heat  exchanger  effectiveness  Is  assumed  only  to  estimate 
an  average  bulk  temperature  for  both  sides  to  determine  the  properties.  Later, 
the  calculated  value  for  the  heat  exchanger  effectiveness  Is  compared  to  the 
Initial  or  previous  value  for  agreement  within  a  specified  tolerance.  After 
determining  the  average  temperatures,  the  fluid  properties,  such  as,  viscosity  u, 
thermal  conductivity  k,  and  specific  heat  at  constant  pressure  Cp  are  ascertained 
which  allows  a  value  for  Prandtl  number  to  be  calculated  (1 .e. ,  NpR  *  Cpu/k). 
Curve-fits  for  the  viscosity  and  thermal  conductivity  of  air  or  a  mixture  of  air 
and  fuel  were  generated  from  tabulated  property  data  published  In  the  open 
literature  [3].  Similar  curve-fits  for  the  specific  heat  of  air  or  a  mixture  of 
air  and  fuel  already  existed  [4].  A  value  for  molecular  weight  of  air  or  mixture 
of  air  and  fuel  Is  required  to  calculate  specific  volume  at  the  core  entrance 
(i.e.,  Vj  »  RyT/pM).  Effects  of  humidity  were  not  Included  In  the  calculations  of 
fluid  properties.  Specific  volumes  are  calculated  using  ratios  of  Inlet  and 
outlet  conditions,  and  then  a  mean  specific  volume  is  determined  for  both  air  and 
gas  sides. 

REYNOLDS  NUMBER 

The  heat  exchanger  flow-stream  mass  velocity  Is  given  by. 


and  expressing  the  Reynolds  number  In  terms  of  G  yields: 
4 r,  G 

it  _  n 
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STANTON  NUMBER,  COLBURN  FACTOR,  AND  FRICTION  FACTOR 

The  friction  factor  f,  and  the  Colburn  factor  j  are  determined  from  the 
tabulated  data  of  reference  1,  corresponding  to  the  selected  fin  surfaces  and  the 
calculated  Reynolds  number.  For  low  Reynolds  numbers  that  are  not  tabulated  In 
reference  1,  the  following  equations  are  used  for  the  f  and  j  factors. 


(14) 


and. 


j 


(15) 


The  value  of  Kj  and  l<2  have  been  determined  by  extending  the  graphs  found  in 
reference  1.  The  Stanton  number  is  extracted  from  the  Colburn  factor  which  Is 
defined  as, 

j  =  N$t  NpR  2/3  .  (16) 

HEAT  TRANSFER  COEFFICIENT 

The  heat  transfer  coefficient  h  is  a  complex  function  of  fluid  properties, 
flow  characteristics,  and  surface  geometries  and  is  defined  by. 


h  "  NST  G  Cp  * 


FIN  EFFECTIVENESS 

The  fin  effectiveness  is  defined  as, 

tanh  m£ 
nf - ml 


where,  m  Is  a  fin  effectiveness  parameter  given  by. 


(17) 


(18) 
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(19) 


m  * 


and  a  Is  the  fin  length  from  root  to  center. 


A 


b 

M7  • 


(20) 


SURFACE  EFFECTIVENESS 

The  total  surface  effectiveness  Is  defined  as, 

no  *  1  '  “T  (1_TV)  * 


(21) 


OVERALL  COEFFICIENT  OF  HEAT  TRANSFER 

Using  the  surface  effectiveness,  heat  transfer  coefficient,  the  total  heat 
transfer  area,  and  the  thermal  wall  resistance;  the  overall  coefficient  of  heat 
transfer  Is  expressed  by 


1  ,  1 
■nr +  TAjK)~njr 

03t  a  g  a  og  g 


i 

k/(a/!2) 


(22) 


NTU  AND  EXCHANGER  EFFECTIVENESS 

The  flow  stream  capacity  rate  mentioned  earlier  is  defined  as, 

C  *  wc  .  (23) 

P 

A  capacity  rate  Is  determined  for  both  sides  of  the  heat  exchanger,  and  a  capacity 
rate  ratio  Is  calculated  for  use  later  in  the  effectiveness  equation  l.e. , 
Equation  (25), 


With  the  aid  of  Equation  (2)  the  number  of  transfer  units  Is  calculated,  which 
with  the  capacity  rate  from  Equation  (23)  are  used  to  calculate  heat  exchanger 
effectiveness.  Recalling  the  counterflow  effectiveness  equation. 
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(25) 


e  * 


.  -NTUU-C  ) 

1  Cre-NTU11-Cr> 


Based  on  the  above  calculated  effectiveness  and  Equation  (la)  and  (lb),  the  air 
and  gas  outlet  temperatures  can  be  calculated. 


PRESSURE  DROP 

The  total  pressure  drop  in  a  compact  heat  exchanger  is  due  to  three  major 
effects:  (1)  air-side  header  pressure  drop,  (2)  core  pressure  drop,  and  (3)  bend 
pressure  drop.  Different  flow  arrangements  can  result  in  lower  or  higher  pressure 
losses.  The  counter-flow  configuration  shown  in  Figure  2  can  result  in  a  uniform 
velocity  across  the  entire  heat  exchanger  core,  [5].  To  keep  the  air-side  header 
pressure  loss  at  a  minimum  and  for  uniform  flow  distribution,  the  dynamic  velocity 
ratio  must  be  kept  a  constant,  [5], 


^  =  0.405  (26) 

q1  7 


where. 


and. 


qi 


V*fe 


(27) 


(28) 


For  the  counter-flow  configuration,  the  total  pressure  drop  die  to  the  headers  on 
the  air-side  is  defined  as 


*ta  q0 

- 1  =  1  -  =  0.595 

qi  q1 


(29) 


To  be  consistent  with  reference  1,  the  header  pressure  drop  must  be  divided  by  the 

Inlet  pressure  pj  , 

3 


n 
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(30) 


By  substituting  the  values  for  q^  and  q0  Into  equation  (26)  the  air-side  velocity 
relationship  can  be  determined. 


V 


o 


.636 


(31) 


Figure  2  -  Heat  Exchanger  Arrangement  Used  for  Pressure  Drop 
and  Weight  Estimations 
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The  core  pressure  drop  equation  for  a  compact  heat  exchanger  Is  complicated 
In  that  it  consists  of  four  components,  [1]:  (1)  entrance  effects,  (2)  flow 

acceleration,  (3)  core  friction,  and  (4)  exit  effects.  The  flow-stream  pressure 
drop  relationship  Is  defined  as, 


*PC  fi2vl 

PT=^VT 


[(K  +l-<£)  +  2(^--  1)  +  f  J-^5. 
c  V1  Ac  V1 


(32) 


entrance  flow  core  exit 

effect  acceleration  friction  effect 


and  for  the  core  friction  parameter. 


(33) 


Definition  of  the  entrance  and  exit  loss  coefficient  are  from  the  literature 
C6].  The  entrance  loss  coefficient,  Kc  Is  given  by. 


K 


c 


<2Kd  - 


(34) 


and  the  exit  loss  coefficient,  Kg  Is  defined  as, 

Ke  *  1  '  2Kd°  +  ^  •  (35) 

Additional  pressure  loss  must  be  considered  since  the  headers  have  matrices.  To 
account  for  the  added  matrix  In  each  header,  the  core  pressure  drop  has  been 
Increased  by  the  percent  of  additional  mean  flow  length  In  the  headers. 
Therefore,  equation  (32)  becomes;  for  the  air-side, 

dp  Lg  Lg 

3  Eq  (32)  x  (— i +  L  +  — j^O/L  (36) 
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and  for  the  gas-side 


Pi 


Eq  (32)  x  (1  + 


(°1  +  °o)  * 
- 2C - J 


(37) 


To  estimate  the  bend  pressure  loss,  it  Is  assumed  that  a  mitre  elbow  is 
similar  to  the  configuration  under  consideration.  The  bend  pressure  loss  is 
calculated  for  the  air-side  and  gas-side  inlet  and  exit  conditions  using  the 
following  expression, 


%  s*  ° 7  2 

T’T** 


(38) 


The  bend  loss  coefficient,  Kjj*  is  curve  fitted  for  the  mitre  elbow  [7]. 
Therefore,  the  total  pressure  drop  of  the  heat  exchanger  for  the  air  side  is  the 
summation  of  equations  (30),  (36),  and  (38)  and  for  the  gas-side  is  the  summation 
of  equations  (37)  and  (38). 

HEADER  DIAMETER  AND  VELOCITY 

Two  Important  parameters  which  deserve  attention  during  the  design  of  heat 
exchangers  are  the  diameters  and  velocities  of  the  air-side  inlet  and  outlet 
headers.  From  equation  (31),  the  velocity  relationship  has  already  been 
determined.  An  Inlet  velocity  must  be  assumed  to  calcuate  an  outlet  velocity.  By 
rearranging  the  equation  of  continuity,  l.e., 

w  =  pAV  (39) 


the  area  of  the  header  pipe  or  opening  Is  determined. 


A  *  -jjjf—  •  (40) 

From  Figure  2,  It  Is  assumed  that  the  area  through  which  the  air  enters  or  exits 
the  manifolds  Is  equal  to  the  area  of  a  circular  pipe  and  Is  given  by. 


A 


(41) 
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Combining  Equations  (40)  and  (41)  the  diameter  of  the  header  is  determined. 


0 


(i 

v-  it 


(42) 


The  conditions  at  which  density  is  defined  depends  on  whether  the  inlet  or  outlet 
header  Is  under  consideration. 

WEIGHT  ESTIMATION 

The  heat  exchanger  weight  estimation  is  comprised  of  five  parts:  (1)  total 
weight  of  air-side  fins,  (2)  total  weight  of  gas  side  fins,  (3)  sum  of  separating 
plates,  (4)  enclosure  weight,  and  (5)  header  weights.  The  calculation  of  fin  and 
plate  weight  is  combined  in  a  single  equation,  but  done  separately  for  the  air- 
side  and  gas-side  because  these  fin  surfaces  may  vary  within  the  regenerator.  The 
basic  equation  for  calculating  fin  and  plate  weight  on  either  side  is  defined  as, 


a 

or 


■  Apm  *  3(1  -  Aft'V/24 


(43) 


A  plate  factor  and  a  fin  factor  ^  in  the  above  equation  are  determined  from 
the  physical  dimensions  of  the  air-side  or  gas-side  fins.  Basically,  the  plate 
factor  ipp,  relates  the  plate's  non-flow  length  to  the  portion  of  non-flow  length 
not  accounted  for  by  fin  thickness.  The  fin  factor,  tyf,  is  calculated  by  dividing 
the  total  area  by  the  extended  fin  area.  Using  the  estimated  fin  and  plate 
weights,  Wa  and  Wg  from  above,  and  core  volume,  V,  an  average  density  for  the  core 
is  determined. 


_  a  g 
pc - V—* 


(44) 


To  help  understand  the  estimated  weight  calculations  of  the  headers  and  enclosure, 
figure  2  is  used.  Note  that  the  headers  are  assumed  to  be  triangular  In  shape. 
The  total  header  volume  which  contributes  weight  to  the  heat  exchanger  is 
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determined  based  on  the  area  of  these  triangle  mil  tipi  led  by  the  core's  non -flow 
length  (l.e.,  square  root  of  frontal  area).  Assuming  header  density  Is  equal  to 
core  density,  header  weight  can  then  be  expressed  as, 

“HD  *  »c  (L9i  Do  L"  *  S  L")/2  •  (45) 

The  enclosure  weight  Is  based  on  the  following  assumptions:  (1)  the  average 
material  density  per  unit  surface  area,  IT  ,  is  approximately  15  lb/ft2  and 
Includes  sheet  metal.  Insulation,  and  supports,  (2)  the  air  manifolds  have  the 
same  perimeter  as  would  the  cl rcumf erences  of  a  circle  with  the  diameter  shown, 
(3)  the  enclosure  Is  comprised  of  a  four-sided  box  with  no  top  and  bottom  and  two 
air  manifolds,  and  (4)  the  headers  form  a  90°  triangle.  From  figure  2  and  the 
preceding  assumptions,  the  enclosure  weight  Is  defined  as, 

UENCL  ■  15  *  (4  L„  L  *  ♦  °o  1  +  L9i  '.,l9  °i  >  •  (46) 

1  0 

Therefore,  the  total  weight  of  the  heat  exchanger  Is  the  summation  of  the  parts 
described  above. 


w  +  w  +  u  +  u 
a  q  HD 


ENCL 


(47) 


PROGRAM  DESCRIPTION 

The  main  program  Is  called  HTER  and  Includes  basic  calculations,  statements 
for  Input  and  output,  and  calls  the  subroutines:  BENDLOS,  STAT,  SURF,  and 
TRANSP.  BENDLOS  calculates  the  bend  loss  coefficient  Kj,*,  used  In  the  pressure 
drop  calculations,  based  on  the  angle  at  which  the  fluid  turns  through.  STAT  Is  a 
data  bank  for  the  friction  factor  and  Colburn  factor  values  from  Kays  and  London 
[1],  and  calls  a  subroutine  INT,  which  Is  an  Interpolation  routine.  SURF  Is  the 
data  bank  for  the  plate  fin  characteristics,  also  from  Kays  and  London  [1], 
(Tables  9-3  a,  b,  c,  and  d).  TRANSP  returns  fluid  properties  of  air  or  combustion 
products.  Including  viscosity,  specific  heat,  and  thermal  conductivity.  Input  to 
the  computer  program  Is  accomplished  by  reading  a  NAMELIST  named  "INDATA".  The 
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NAMELIST  INDATA  Identifies  a  succeeding  list  of  input  variables  which  can  be  Input 
without  speclflng  format.  Included  In  the  program  is  a  BLOCK  DATA  DEFAULT  for  the 
NAMELIST  INDATA,  which  defines  a  default  value  for  each  input  variable.  More 
about  the  NAMELIST  statement  and  default  values  Is  given  In  the  section  titled 
INPUTS.  The  main  program  and  subroutines  are  written  In  FORTRAN  IV  language  for 
the  CDC  6600  or  6700  computers. 

The  computer  program  has  the  following  limitations  and  features,  [2]: 

1.  The  program  assumes  even  distribute  >n  of  flow  across  the  flow  cross 

section;  i.e.,  the  mass  flow  rate  through  any  channel  Is  simply  the  mass  flow 
for  sectional  area  of  the  channel.  If  this  assumption  Is  not  true, 

effectiveness  drops  rapidly.  Even  flow  distribution  into  the  heat  exchanger 
is  not  easy  to  achieve,  especially  In  a  short  flow  length,  low  pressure  drop 
exchanger  [8]. 

2.  Effects  of  fouling  and  fin  deformation  are  not  Included  in  the 
program. 

3.  An  additional  heat  transfer  (~2  to  5%)  in  a  counterflow  exchanger  is 
obtained  from  the  cross-counterflow  In  the  headers.  This  was  not  considered 
In  the  program.  In  a  pancake-shaped  counterflow  exchanger,  the  error  in  this 
assumption  becomes  appreciable. 

4.  All  plate-fin  surfaces  given  in  Figures  9-3  to  9-7  of  Kays  and 
London  [1]  are  Included  in  the  program  except  for  4.0  fins/in  (plain). 

The  following  sections  describe  the  necessary  input  variables  and  resultant 
output  variables  of  the  computer  program.  Several  sample  cases  of  Input  data  are 
given  Illustrating  format,  as  well  as,  their  corresponding  output.  Also  described 
are  error  messages  that  may  occur  during  the  execution  of  the  program.  A  complete 
listing  of  the  computer  program  Is  Included  In  the  report  as  APPENDIX  A,  and  the 
main  program  variables  are  defined  In  APPENDIX  B. 
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INPUTS 

Physical  heat  exchanger  core  parameters  (l.e.,  lengths,  fin  types,  frontal 
areas,  etc.),  air-side  Inlet  header  velocity,  and  heat  exchanger  flow  conditions 
(l.e.,  pressure,  temperatures,  mass  flows  and  fuel-air  ratio),  are  provided  to  the 
program  In  the  form  of  a  NAMELIST  called  INDATA.  The  Input  conditions  are 
required  to  be  in  U.S.  units.  A  set  of  default  values  have  been  Included  in  the 
program  and  Include  engine  cycle  data  from  figure  10  of  reference  9.  NAMELIST 
INDATA  variables  are  as  follows  (with  the  default  values  shown  to  the  right): 

TYPA  -  air  side  fin  type  Default  Value  =  1 

*  1  -  plain  fin 

3  2  -  louvered  fin 
■  3  -  strip/offset  fin 

*  4  -  wavy  fin 

NSA  -  air  side  surface  number  *  7 

*  1  to  18  for  plain  fin 

*  1  to  14  for  louvered  fin 

*  1  to  12  for  strip/offset  fin 

*  1  to  3  for  wavy  fin 

TYPG  -  gas  side  fin  type  3  1 

=  1  -  plain  fin 
=  2  -  louvered  fin 
3  3  -  strip/offset  fin 
3  4  -  wavy  fin 

NSG  -  gas  side  surface  number  3  7 

3  1  to  18  for  plain  fin 
3  1  to  14  for  louvered  fin 
3  1  to  12  for  strip/offset  fin 
3  1  to  3  for  wavy  fin 

RLENI  -  Initial  flow  length,  ft  3  3.0 
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RLI 

-  Incremental  flow  length,  ft 

= 

1.0 

NL 

-  number  of  length  Iterations 

= 

5 

IAFRA 

-  Initial  frontal  area,  ft2 

* 

25.0 

AFRAI 

-  value  of  frontal  area  Increase  per  Iteration,  ft2 

= 

25.0 

NA 

-  number  of  frontal  area  Iterations 

= 

4 

WA 

-  air  side  mass  flow  rate,  lbm/sec 

= 

90.0 

PINA 

-  air  side  inlet  pressure,  lb^/ln2 

= 

116.4 

TINA 

-  air  side  inlet  temperature,  °R 

= 

1040.5 

PEXG 

-  gas  side  outlet  pressure,  1bf/in2 

= 

14.9 

WG 

-  gas  side  mass  flow  rate,  lbm/sec 

s 

101.45 

TING 

-  gas  side  Inlet  temperature,  °R 

— 

1646.4 

FAR 

-  fuel-to-alr  ratio 

= 

0.0145 

V  INPUT* 

-  header  velocity,  ft/sec 

s 

90.0 

SAMPLE  CASES 

Included  In  this  manual  are  eight  (8)  sample  cases  to  Illustrate  using 
NAMELIST  INDATA  to  define  sets  of  Input  data.  The  data  cards  for  these  sample 
cases  are  shown  In  Figure  3.  INDATA  variable  sets  must  begin  with  a  SINDATA 


*  The  user  Is  advised  to  keep  Inlet  air  velocity  below  100  ft/sec  to  ensure  a 
design  with  good  flow  distribution  and  low  losses 


19 


PAS  82-41 


(where  the  $  Is  located  In  the  second  column)  and  each  set  ends  with  a  $.  The 
first  data  card  has  the  IINDATA  $,  therefore,  the  default  values  are  used  to 
calculate  the  results  shown  as  TABLE  C-l  In  APPENDIX  C.  A  set  of  data  Items  may 
consist  of  any  subset  of  the  variables  names  In  NAMELIST  INDATA.  The  value  of 
variables  not  Included  In  the  subset  on  a  following  card  remains  unchanged. 


z 


/  jINOATA  TYPA»3 ,NSA«10,TYPG»3 ,NSG»1 ,RLSNI»1.0  % 
f  SI NO AT A  NSA»5,NSG»2,RLENI»14.Q,ML»8,NA«3  $ 

/  SINOATA  TYPA«l,NSA»ll,NSQ7fRLENI»2.0  % _ 

/  SINOATA  NS6»5,RLENI»4.Q  $ 


IINDATA  TYPA-2  % 


/ 


/  SINOATA  PINA»223.6,TINA»878.9,WG»101.75  JING-142S.8,FAR»oTq17S  S 
/  SINOATA  PINA-111. 9,TINA»730.5,WG»101.49,TING»1641.3,FAR*0.0149  S 
/SINOATA  %  ™ 


Figure  3  -  Sample  Cases  For  Namelist  INDATA 

The  output  consists  of  title  headings.  Input  variables,  and  computed 
results.  APPENDIX  C  contains  the  correspondl ng  output  for  the  sample  cases  shown 
In  Figure  3.  Each  data  card  produced  a  page  of  output.  The  first  line  of  the 
output  Is  a  title  heading:  "CORE  HEAT  TRANSFER  SURFACES".  The  next  six  lines  of 
output  give  the  heat  transfer  surface  characteristics  (l.e. ,  hydraulic  radius, 
compactness,  plate  spacing,  etc.).  "HEAT  EXCHANGER  CONDITIONS"  Is  the  next  title 
heading  and  given  on  the  following  four  lines  are  the  mass  flows,  pressures, 
temperatures,  and  fuel -air  ratio  for  the  air  and  gas  streams.  These  values  are 
the  INDATA  variables.  The  next  title  heading  Is  "HEADER  DESIGN  DETAILS,"  and 
given  on  the  next  four  lines  are  the  air-side  Inlet  diameter  and  velocity.  The 
computed  values  for  the  remaining  output  variables  at  various  combinations  of  core 


20 


PAS  82-41 


f 


counterflow  lengths  and  frontal  areas  comprise  the  rest  of  the  output.  With  each 
new  frontal  area  Iteration,  new  column  descriptors  are  printed.  Contained  In  the 
tabulated  values  are  as  follows:  counterflow  length;  heat  exchanger  core  volume; 
frontal  area;  air  side  outlet  pressure;  gas  side  Inlet  pressure;  air,  gas,  and 
total  pressure  drops  In  percent;  air  and  gas  side  outlet  temperatures;  overall 
enclosure  height;  number  of  transfer  units;  heat  exchanger  effectiveness  In 
percent;  estimated  weight;  air-side  exit  header  diameter  and  exit  header  velocity. 

The  results  shown  in  Table  C-l  (see  APPENDIX  C)  were  generated  assuming  plain 
fins,  11  fins  per  Inch,  and  0.25  inch  plate  spacing  on  both  the  air  and  gas 
sides.  Counter-flow  lengths  are  varied  from  3  to  7  feet  on  1  foot  Increments  and 
total  frontal  flow  area  Is  varied  from  25  to  100  ft2  in  25  ft2  increments.  Air 
and  gas  Inlet  conditions  are  based  on  the  gas  turbine  cycle  data  shown  in  Figure 
10  in  reference  9.  An  air-side  Inlet  header  velocity  of  90  ft/sec  Is  assumed  in 
all  sample  cases.  Although  the  physical  characteristics  of  the  heat  exchangers 
remain  the  same,  air  and  gas  side  Inlet  conditions  were  changed  in  Tables  C-2  and 
C-3  to  reflect  gas  turbine  cycle  conditions  shown  as  Figures  11  and  12, 
respectively,  in  reference  9.  Air  and  gas  side  inlet  conditions  remain  the  same 
for  Tables  C-3  through  C-8  while  physical  characteristics  of  the  heat  exchangers 
are  varied.  Table  C-4  shows  the  effect  of  changing  the  air-side  fin  type,  from 
plain  to  louvered,  but  same  fins/in,  has  on  the  performance  of  the  heat 
exchanger.  Table  C-5  varies  the  gas-side  flns/ln  from  11.1  to  6.2  and  changes  the 

Initial  flow  length  from  3  feet  to  4  feet.  In  Table  C-6  variations  include  the 

following:  (1)  air-side  fin  type  from  louvered  to  plain,  (2)  air-side  and  gas-side 
flns/ln  from  11.1  to  19.86  and  6.2  to  11.1,  respectively,  and  (3)  redefines  the 
Initial  flow  length  from  4  feet  to  2  feet.  Table  C-7  varies  the  air-side  and  gas- 
side  flns/ln  from  19.86  to  6.2  and  11.1  to  3.01,  respectively,  and  varies  the 

Initial  flow  length  from  2  feet  to  14  feet.  Also,  varied  in  Table  C-7  are  the 

number  of  length  Iterations  and  frontal  area  Iterations  from  5  to  8  and  4  to  3, 
respectively.  Table  C-8  varies  the  following  assumptions:  (1)  fin  types  for  both 
sides  (air  and  gas)  from  plain  to  strip/offset,  (2)  the  Initial  flow  length  from 
14  feet  to  1  foot,  and  (3)  flns/ln  on  the  air-side  and  gas-side  from  6.2  to  19.82 
and  3.01  to  11.1,  respectively. 
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ERROR  MESSAGES 

There  are  two  basic  error  messages  that  may  occur  in  the  execution  of  the 
program.  These  messages  are  presented  In  this  section  along  with  reasons  for 
their  occu ranee. 

(1)  REYNOLDS  NUMBER  OUT  OF  RANGE  OF  PROGRAMMED  TABLES  -  The  Reynolds  number 
on  the  air  or  gas  side  being  too  large  causes  this  error  message  to  occur. 
The  reason  why  it  occurs  is  that  only  the  values  from  reference  1  were 
tabulated  in  the  computer  program  and  have  an  upper  limit.  The  lower  limit 
has  been  already  discussed  in  the  Method  of  Analysis. 

(2)  TRANSP  INPUT  OUT  OF  RANGE  -  The  reason  for  this  error  message  is  that  the 
entering  fuel-air  ratio  to  the  subroutine  TRANSP  is  too  large  or  too  small.  A 
maximum  value  of  0.034826,  and  a  minimum  value  of  less  than  0.0  are  associated 
with  the  curve-fits  for  certain  thermodynamic  properties.  Also,  if  the 
temperature  is  less  than  500  °R  or  greater  than  2000  °R  the  error  message  will 
occur.  Results  obtained  within  these  ranges  agree  with  the  data  from 
reference  3  within  ±  1*. 
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♦DECK  HTER 

PROGRAM  HTER ( I NPUT , OUTPUT , TAPES- I NPUT , T APE6-0UTPUT ) 

C 

C****#**##**##***#**********##****#*4Ht*********#*******#*****#**#«***** 


HEAT  EXCHANGER  TYPE  TO  BE  SPECIFIED  AS  FOLLQWS- 
WHERE  TYPE: 

TYPE  1  -  PLAIN  FIN 
TYPE  2  -  LOUVERED  FIN 
TYPE  3  -  STRIP/OFFSET  FIN 
TYPE  4  -  WAVY  FIN 


SURFACE  NUMBER 
SURFACE  NUMBER 
SURFACE  NUMBER 
SURFACE  NUMBER 


TO  IB 
TO  14 
TO  12 
TO  3 


C******#*************##**#***##*#****#**************#***********'******* 

c 

INTEGER  TYPA,  TYPG 
C 

REAL  I LA,  ILG,  K,  KA,  KBINA,  KBEXA,  KBING,  KBEXG,  KG,  KCA,  KCG. 

1  KDA,  KDG,  KEA,  KEG,  LA,  LG,  LHGING,  LHGEXG,  MA,  MG,  MAL.  MGL, 

2  MUA,  MUG,  NARA,  NARG,  NE,  NPRA,  NPRG,  NRA,  NRG,  NSTA,  NSTG, 

3  NTU,  XAfiA 
C 

COMMON/ I NPUT / TYPA , NSA , TYPG , NS6, RLENI , RLI , NL, I AFRA , AFRA I , NA , WA , P I NA 
2  , TINA. PEXG.WG, TING, FAR. VINPUT 

C 

NAMELIST/ INDATA /TYPA, NSA, TYPG, NSG, RLENI , RLI , NL, I AFRA, AFRA I , NA, WA, 

2  PINA, TINA, PEXG, WG, TING. FAR, VINPUT 

C 

DATA  RHO,  PI,  RU/  4B5. ,  3.1416.  640.1  / 

DATA  A, K, GC, ERROR/ . 01 , 12. ,32.2, .001/ 

DATA  E, DELPA, DELPG  /.5,  0.01,  0.03/ 

C 

33  READ (5, INDATA) 

IF (EOF (5))  999,65 
C 


C#**#*******###*  HEAT  EXCHANGER  CORE  DIMENSIONS  ********************* 

C 

63  CALL  SURF (TYPA, NSA, AXA,BXA,SFA, BA, RHA.DELA, BETA, FRA, WFA,WPA) 
CALL  SURF (TYPG, NSG, AXG, BXG, SFG, BG, RHG, DELG, BETG, FRG, WFG, WPG) 


»# _ II 

II  _  II 


VINA-VINPUT 

RHO I N A-P I NA / T I N A / RU* 1 728 . 

DINA— SORT (4. *WA/RHOINA/VINA/PI ) 

WRITE (6, 110)  TYPA, AXA.NXA.SFA, TYPG, AXG, NXG, SFG, BA, BG.RHA, RHG, 

DELA, DELG, BETA, BETG, FRA, FRG, WA, WG, PINA. PEXG, TINA,  TING 
, FAR, FAR 

( 1H1 , "  CORE  HEAT  TRANSFER  SURFACE" , 18X , "AIR-SIDE" , 29X , 
"GAS-SIDE"  //  20X,  "TYPE  AND  FIN  DETAIL",  3X.I1, 

F6.4,  II,  F5.2,20X,  II,  F6.4,  II, 

,  F5.2  /  20X ,  "PLATE  SPACING",  12X,  F6.4,  2X,  "IN",26X, 
F6.4,  2X ,  "IN"  /20X ,  "HYDRAULIC  RADIUS",  SX,  F7.5,  "  FT", 
25X , F7. 5,  "  FT"  /20X,  "FIN  THICKNESS",  12X,  F6.4,  2X, 

" IN" , 26X ,  F6.4,  2X ,  "IN"  /20X,  "COMPACTNESS",  13X,  F7.1, 

2X , "SOFT /CUFT" , 15X.F10. 1,  2X,  "SQFT/CUFT"  /  20X,  "FIN/TOTAL", 
"  AREA",  11X,  F6.4.  2X,  "FT/FT", 23X,  F6.4,  2X,  "FT/FT"// 

"  HEAT  EXCHANGER  CONDITIONS" , 1BX, "AIR-SIDE  INLET", 24X, 
"GAS-SIDE  INLET", 10X, "GAS-SIDE  EXIT"//20X, "MASSFLOW"  ,  16X, 
F6.2, "  LB/SEC", 23X.F6. 2, "  LB/SEC", 17X, "-"/20X, "PRESSURE" 
, 16X , F6. 2, "  PSIA" , 32X , "-" , 1SX , F6. 2, "  PSIA"/20X, "TEMPERATURE", 
12X, F7. 2, "  DEG  R" , 23X , F7. 2, "  DEG  R" . 18X, "-"/20X, "FUEL-" 


2 
3 

110  FORMAT 
2 

3 

4 

5 

6 

7 

8 
9 


5"AIR  RATIO" , 12X , "0. 0" , 34X , F6. 4, 20X , F6. 4//"  HEADER  DESIGN  DETAILS" 
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6  /) 

WRITE (6, 115)  VINA, DINA 

115  FORMAT (20X, "INLET  AIR  HEADER  DIAMETER  SIZED  FOR  INLET", 

2  "  AIR  VELOCITY  =",F8.2,"  FT/SEC" /20X , "OUTLET  ", 

3  "AIR  HEADER  DIAMETER  SIZED  FOR  UNIFORM  FLOW  DISTRIBUTION  " 

4  "AND  MINIMUM  HEADER  LOSS" /20X ," INLET  DIAMETER  =  ", 

5  F6.2,"  FT"/20X, "EXIT  AIR  DIAMETER  AND  VELOCITY  GIVEN  " 

6  "BELOW"/) 

DO  700  1=1, NA 
11=1-1 

AFRA- 1 AFRA+ 1 1 *AFRA I 
WRITE (6, 130) 

130  FORMAT <1H0,"  LENGTH  VOLUME  AREA  P-A-EX  P-G-IN  DPA  DPG 
1  DPT  T-A-EX  T-G-EX  HEIGHT  NTU  EFFECT  WEIGHT  HEADE 
2R  VELOCITY"/"  FT  CU  FT  SOFT  PS I A  PS I A  PCT 

3PCT  PCT  DEG  R  DEG  R  FT  ",10X,"PCT  LBS  DIA 

4FT  FT/SEC"/) 

DO  600  L-1,NL 
Ll-L-1 

RLEN=RLENI+RLI*L1 
VOL=AFRA*RLEN 
PEXA=PINA*  < 1-DELPA) 

P I NG=PEXG*  < 1 +DELPG ) 

C 

C***************  HEAT  TRANSFER  AND  FREE  FLOW  AREAS  ****************** 

C 

PAA=BA/ (BA+BG+2. *A) 

PAG=BG/ (BA+BG+2. *A) 

ALHA=BETA*PAA 

ALHG=BETG*PAG 

AFRG=AFRA 

AA=ALHA*VOL 

AG=ALHG*VOL 

S I GA= ALH A*RHA 

SIGG«ALHG*RHG 

AC A=S I GA*AFRA 

ACG=S I GG* AFRG 

FLA-VOL/AFRA 

XNCFL-SQRT(AFRA) 

NXA-BXA 
NXG-BXG 
CMIN=1 . 

CA-1. 

CG-1. 

C 

C***************  AVERAGE  CORE  FLUID  PROPERTIES  ************************* 
C 

150  TEXA-  E* (TING-TINA) *CMIN/CA+TINA 
JAVA- (TEXA+TINA) *. 5 

CALL  TRANSP  (TAVA,  0.,  CPA,  KA,  MU A,  MA) 

TEXG-  E* (TINA-TING) *CMIN/CG+TING 
TAVG- (TEXG+TING) *.  5 

CALL  TRANSP  (TAVG,  FAR.  CPG,  KG,  MUG,  MG) 

C 

C**<mhhhh»*******  AVERAGE  CORE  REYNOLDS  NUMBERS  *####*##*#**•*##*■*##*#***# 
C 

GA-WA/ACA 
NRA=4. *RHA*GA/MUA 
GG-WG/ACG 
NRG-4. #RHG*GG/MUG 
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C***************  STANTON  NUMBER, COLBURN  FACTOR, AND  FRICTION  FACTOR  **** 
C 

CALL  STAT  (TYPA,  NSA,  NR A,  COLBFA, FA) 

IF  (FA.  EQ.  0.  )  GO  TO  700 
NPRA»CPA*MUA/KA 
NSTA«COLBFA/NPRA**.  666 
HA*NST A*GA*CPA*3600 . 

CALL  STAT  (TYPQ,  NSG,  NRG,  COLBFG,  FG) 

IF  <FG. EQ. 0. )  GO  TO  700 
NPRG”CPG*MUG/KG 
NSTG=COLBFG/NPRG**. 666 
HG”NSTG*GG*CPG*3600. 

C 

C***************  FIN  EFFECTIVENESS  ****************************** 

C 

MAL=SORT<  <2. *HA) / <K*DELA/12. ) ) 

MGL=SQRT  (  (2.  *HG)  /  OODELG/12.  )  ) 

LA= (BA/12. ) /2. 

LG=(BG/12. ) /2. 

MAL“MAL*LA 
MGL=MGL*LG 
ETAFA-TANH(MAL) /MAL 
ETAFG=TANH (MGL) /MGL 
C 

C***************  SURFACE  EFFECTIVENESS  ********************** 

C 

ETAOA-1 . -FRA* ( 1 . -ETAFA) 

ETAOG”l . -FRG* < 1 . -ETAFG) 

C 

C***************  OVERALL  COEFFICIENT  OF  HEAT  TRANSFER  ************* 

C 

RA-1. / <ETAOA*HA) +1 . / ( < AG/AA) *ETAOG*HG) +1 . / <K/ (A/12. ) ) 

UA=1. /RA 
C 

C***************  INLET  AND  EXIT  LOSS  COEFFICIENTS  ******************** 

C 

CCA” . 6 1 00000000 1 - . 1 444294507 1 *S I GA+ 1 . 0080347435*S I GA**2 
CCA=CCA-1 . 73 1 7560083*S I GA**3+ 1 . 1559407939*SIGA**4 
CCG” . 6 1 00000000 1 - . 1 444294507 1 *S I GG+ 1 . 00B0347435*S I GG**2 
CCG-CCG-1 . 73 1 75600B3*S I GG**3+ 1 . 1559407939*SIGG**4 
NARA=NRA*l.E-4 
NARG=NRG*1 . E-4 

KDA=1. 10639601 04-. 13322445533*NARA+. 1 188542B625*NARA**2 
KDA-KDA-. 033170530592*NARA**3 

KDG-1. 10639601 04-. 13322445533*NARG+. 1 1885428625*NARG**2 
KDG=KDG- . 033 1 70530592*NARG**3 
KCA” <1.-2. *CCA+CCA**2* <2. *KDA-1 . ) ) /CCA**2 
KCG- ( 1 . —2. *CCG+CCG**2* (2. *KDG-1 . ) ) /CCG**2 
KEA-  1.-2. *KDA*S I G A+S I GA**2 
KEG”  1.-2. *KDG*S I GG+S I GG**2 
C 

C***************  PRESSURE  DROPS  ****************************** 

C 

RH0EXA-PEXA/TEXA/RU*1728. 

RHO I NG-P I NG /T I NG /RU* 1 728 . 

RHOEXG»PFXG/TEXG/RU*1728. 

VEXA-. 636*VINA*SQRT (RHOINA/RHOEXA) 

DEXA*SQRT  <4. *WA/RHOEXA/VEXA/PI ) 

HINA-RHOINA/2. /GC*VINA**2 
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DELPAH- . 595*H I NA/P INA/ 1 44. 

I LA* 1 . -S I GA**2+KCA 

ILG-1 . -SIGG**2+KCG 

SPVA=PINA/PEXA*TEXA/TINA 

SPVG=PING/PEXG*TEXG/TING 

ELA= ( 1 . — S I GA**2— KEA ) *SPVA 

ELG= < 1 . -SIGG**2-KEG)*SPVG 

SPVAM-2 . *P I NA/ <  P I NA+PEX A ) *TAVA/T I NA 

SPVGM—2 . *P I NG / (PING+PEXG) *TAVG/TING 

cfa=fa*aa/aca*spvam 

CFG»FG*AG/ACG*SPVGM 
FAA=2. * (SPVA-1 . ) 

FAG-2. * (SPVG-1 . ) 

TLA-ILA+FAA+CFA-ELA 
TLG- I LG+FA6+CFG-ELG 
LHG I  NG-SQRT  ( XNCFL**2-DEX A**2 ) 

LHGEXG-SQRT ( XNCFL**2-DINA**2) 

ANGEXG-ATAN (DINA/LHGEXG) 

ANGING-ATAN  (DEXA/LHGING) 

ANGINA-PI/2. -ANGEXG 
ANGEX A-P 1/2. -ANG I NG 
HFXG-1 .  +  < (DINA+DEXA) /2. /FLA) 

HFXA- (LHGING/2. +FLA+LHGEXG/2. ) /FLA 

DELPACXGA/144.  /PINA)  **2/2.  /GC*1545.  /MA*TINA*TLA*HFXA 

DELPGC— (GG/ 144. /PING)**2/2. /GC*1545. /MG*T I NG*TLG*HFXG 

AHCM I NA=S I GA*D I NA*XNCFL 

AHCM I  NG-S  I  GG*XNCFL*LHG  I  NG 

V I NAH 1 -WA/RHO I N A / AHCM I N A 

VINGH1-WG/RH0ING/AHCMING 

VINAC1-VINAH1*C0S( ANGINA) 

V I NGC 1 = V I NGH 1 *COS  <  ANG I NG ) 

CALL  BENDLOS  (ANGINA, KB INA) 

CALL  BENDLOS  ( ANGING, KBING) 

VINAM-SQRT  <  (VINAH1**2+VINAC1**2> /2. > 

VINGM-SCRT  <  < VINGH1**2+VINGC1**2) /2.  ) 

DELPAB 1 =RHO I NA*KB INA/2. /GC*V I NAM**2 

DELPGB 1 -RHO I NG*KB ING/2. /GC*V I NGM**2 

AHCMEXA-SIGA*DEXA*XNCFL 

AHCME  XG-S I GG*LHGE  X  G*  X  NCFL 

VE  X  AH2-WA /RHOE  X  A / AHCME  X  A 

VEXGH2-WG/RH0EXG/AHCMEXG 

VE  X  AC2- VE  X  AH2*C0S ( ANGE  X  A ) 

VEXGC2-VEX0H2*COS (ANGEXG) 

CALL  BENDLOS  ( ANGEXA, KBEXA) 

CALL  BENDLOS  (ANGEXG, KBEXG) 

VEXAM-SQRT( (VEXAH2**2+VEXAC2**2) /2.  ) 

VEXGM-SQRT( <VEXGH2**2+VEXGC2**2> /2.  > 

DELPAB2»RH0EXA*KBEXA/2. /GC*VEXAM**2 
DELPGB2»RH0EXG*KBEXG/2. /GC*VEXGM**2 
DELPAB- ( DELPAB 1 +DELPAB2) /P I NA/ 1 44 . 

DELPGB- ( DELPGB 1 +DELPGB2 ) /P ING/ l 44 . 

DELPA-DELPAC+DELPAH+DELPAB 
DELPG-DELPGC+DELPGB 
PEXA-PINA* ( 1 . -DELPA) 

PING— PEXG* ( 1 . +DELPG) 

PCDELPA- 1 00 . *DELPA 
PCDELPG- 1 00 . *DELPG 

C 

C***************  NTU  AND  HEAT  EXCHANGER  EFFECTIVENESS  ***************** 

C 
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CA»WA*CPA*3600. 

CG-WG*CPG*3600. 

CMIN-AMIN1 (CA.CG) 

CR-AMIN1 (CMIN/CA,  CMIN/CG) 

NTU«AA*UA/CMIN 

C  COUNTER-FLOW  EFFECT I VENESS-NTU  RELATIONSHIP 
X»EXP(-NTU*(1.-CR) ) 

NE-  < 1 . — X) / (l.-CR*X) 

IF  (ABS  (NE/E-1.)  .LT.  ERROR)  GO  TO  300 

E-NE 

GO  TO  150 

300  TEXA-NE* (TING-TINA) *CMIN/CA+T1NA 
TEXG-NE* (TINA-TING) *CMIN/CG+TING 
DELPT  -DELPA+DELPG 
PCDELPT* 1 00 . * DELPT 
PCNE* 1 00 . *NE 
C 

C***************  WEIGHT  CALCULATIONS  OF  THE  HEAT  EXCHANGER  *********** 

C 

WTA-AA*RHO* (FRA*DELA*WFA  +  A*  ( 1 .  -FRA)  *WPA>  /24. 
WTG=AG*RHO*(FRG*DELG*WFG  +  A* ( 1 . -FRG) *WPG> /24. 

WPLA- 1 5 . * ( 4 .  *XNCFL*FLA+P I *XNCFL* ( D I NA+DEX A  > +LHG I NG*DEX A 
2  +LHGEXG*DINA> 

WIE= (WTA+WTG) /VOL* (LHGING*DEXA*XNCFL+LHGEXG*DINA*XNCFL) /2. 

WHXT—WT  A+WTG+WPLA+W I E 
OVALHT-D I NA+DEX A+FLA 

WR I TE  < 6 . 500 ) RLEN , VOL , AFRA , PEXA, PING, PCDELPA , PCDELPG , PCDELPT , TE  X A , 

2  TEXG, OVALHT, NTU, PCNE, WHXT, DEXA, VEXA 

500  FORMAT (F7. 1 , F9. 1 , FB. 1 , 2F9. 2, 3F6. 2, 2F9. 1 , F9. 2. FB. 2, F8. 2, F9. 1.F9.2 
2  ,  FI  1 . 2) 

VINA-VINPUT 
600  CONTINUE 
700  CONTINUE 
GO  TO  35 
999  STOP 
END 

*DECK  STATIS 

SUBROUTINE  ST AT  (TYPE,  NN,  RE,  NST ,  F) 

#*#**###*##*##***#*#*####**#*#####*###*###*4*****#***#*####*##* 

SUBROUTINE  STAT  RETURNS  STANTON  NUMBERS  AND  FRICTION  FLOW  DATA 
FOR  THE  TYPE  HEAT  EXCHANGER  SPECIFIED 
**##**************#********************************#*******+*** 
INTEGER  TYPE 

REAL  NR, NST, I R , J R , KR , LR , IS, JS, KS, LS, MS, NS, MR 

DIMENSION  AR<4, 18) ,BR<4, 18) ,CR<4, 18) ,DR(4, 18) ,ER(4, IB) ,FR(4, 18) , 
*IR(4, 18) , JR (4, IB) ,KR<4, 18) ,LR<4, 18) ,MR<4, 18) ,NR<4, 18) ,0R(4, 18) , 

*PR <4, 18) ,QR(4, 18) ,SR<4, 18) ,TR(4, 18) ,UR<4, IS) ,VR<4, 18)  ,WR(4,  18)  , 

*XR (4, 18) ,YR<4, 18) ,AS<4, 18) ,BS(4, 18) ,CS<4, 18) ,DS<4, 18) ,ES(4, 18) , 

#FS  <4, 18) ,GS(4, 18) ,HS(4, 18) , IS (4, 18) , JS<4, 18) ,KS<4, IS) ,LS(4, 18) , 

*MS (4, 18) .NS (4, 18) ,0S(4, 18) .PS (4, 18) ,QS<4, 18) ,RS<4, 18) ,TS(4, 18) , 

*US (4, 18) ,VS(4, 18) ,WS(4, 18) ,YS<4, 18) ,XS<4, 18) ,ZS(4, 18) 

*, S ( 18) ,S1 (18) , S2 ( 18) , S3 (18) ,C(14) , D ( 14) 

DATA ( (S ( I ) , 1-1, 18)- 

*2.0,3.01,3.97.5.3,6.2,9.03, 11. 1. 11. 11, 14.77, 15.08. 19.86. 10.27, 
*11.94,12.00,16.96,25.79,30.33,46.45) 

DATA ( ( <AR(N, J) , J-l, 18)  ,N-1,4)- 

*60000 .  ,  50  *'  1 .  ,  40000 .  ,  30000 .  .  25000 .  ,  20000 .  ,  1 5000 .  .  1 2000 .  ,  1 0000 .  ,  800 
*0. ,6000. ,w  100. ,4000. ,5*0. , .00228, .00237, .00248, .00264, .00274, .0028 
*8 , . 00305 , . 00320 , . 00333 , . 00347 , . 00363 , . 00373 , . 00379 , 5*0 . , . 00549 , . 00 
*562, .00579, .00601, .00616, .00638, .00672, .00703, . 00734, . 00778, .00847 
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*, .00904, .01023,5*0. ,0.0591,25.00, 16*0.0) 

DATA  < <  <BR(N, J) , J  =  l, 18) ,N-1,4)» 

*45000.  ,40000.  .30000.  ,25000.  .20000. , 15000. , 12000. , 10000. ,0000. .6000 
*. ,5000. ,4000. ,3000. ,5*0. , .00233, .00239, .00254, .00264, .00277. .00295 
*, .00310, .00322, .00336, .00355, .00366, .00373, .00368,5*0. . .00602, .006 
*08, . 00630, . 00645, . 00667, . 00700, . 00732, . 00762, . 00808, . 00886, . 00950. 
*.01045, .01190,5*0. ,0.8591,25.00, 16*0.0) 

DATA ( ( ( CR ( N ,  J  ) ,  J  =  1 , 18) , N- 1,4)= 

*35000. , 30000. , 25000. , 20000. , 15000. , 12000. , 10000. , 8000. , 6000. . 5000. 
*,4000. ,3000. ,2500. ,5*0. , .00246, .00254, .00263, .00276, .00291, .00302. 
*.00316, .00330, .00348, .00357, .00367, .00367, .00357,5*0. , .00595. .0060 
*5, .00620, .00638, .00667, .00695, .00720, .00761, .00826, .00880, .00963, 
*.01110, .01230,5*0. ,0.8591,25.00, 16*0.0) 

DATA  <  (  (DR(N, J) , J  =  l, 18) ,N=1.4)= 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. .2500. .2000. , 1500. , 1200. , 1000. 
* ,  800 .  ,  600 .  ,  5*0 .  , . 00373 ,. 00397 , . 00427 , . 00448 , . 00477 , . 005 15,. 00535 . . 
*00554, . 00571 , . 00606, . 00654, . 00728, . 00851 , 5*0. , . 00764, . 00806, . 00870 
*, .00913, .00978, .0108, .0115, .0127, .0146, .0167, .0189, .0228, .0299,5*0 
*. ,0.7839, 18.36, 16*0.0) 

DATA  < <  <ER(N, J) , J-l, 18)  ,N«1,4>* 

*12000.  , 10000.  ,  8000.  ,  6000. , 5000. , 4000. , 3000. , 2500. , 2000. , 1500. . 1200 
*. , 1000. ,800. ,5*0. , .00303, .00310, .00317, .00325, .00330, .00333, .00326 
*, .00301, .00312, .00371, .00435, .00496, .00581,5*0. , .00708, .00735, .007 
*68, .00807, .00838, .00875, .00923, .00958, .0103, .0127, .0152, .0176, .021 
*1 , 5*0. ,0.6251, 17.24, 16*0. 0) 

DATA ( <  (FROM,  J)  ,  J  =  l,  18)  ,N=1,4>  = 

*15000. , 12000. , 10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 150 
*0. , 1200. , 1000. ,800. ,4*0. , .00255, .00265, .00273, .00283, .00296. .00304 
*, .00310, .00310, .00318, .00347, .00421, .00499, .00575, .00692,4*0. , .007 
*08, .00740, .00763, .00799, .00842. .00870, .00903, .00980, .0106, .0122, 
*.0152, .0182, .0214, .0262.4*0. ,0.7233,21.40, 16*0.0) 

DATA  ( ( (IR(N, J) , J-l, 18) ,N=1,4>* 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. , 1000. 
*,800. ,600. ,500. ,4*0. , .00314, .00333, .00356, .00372. .00390, .00412, .00 
*424,  .00436,  .00444, .00471, .00515, .00599, .00733, .00840,4*0. . .00878, . 
*00923, .00971, .00991, .0103, .0112, .0119, .0139, .0149, .0169, .0190, .022 
*8, .0294, .0350,4*0. ,0.6471, 17.80, 16*0.0) 

DAT  A ( < ( J  R  ( N ,  J  ) ,  J  =  1 , 18) , N= 1 , 4  >  = 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. , 1000. 
*,800. ,600. ,500. ,4*0. , .00288, .00303, .00324, .00338, .00353, .00368, .00 
*373 ,  .  00375 ,  .  00420 , . 00505 , . 00586 , . 00704 , . 00890 ,.0103, 4*0 . , . 00763 , . 0 
*0807, .00862, .00900, .00958, .0105, .0112, .0119, .0137, .0166,  .0198.  .024 
*3, .0319, .0380,4*0. ,0.7800, 19.19, 16*0.0) 

DATAC ( <KR(N. J) , J-l, 18) ,N»1,4) - 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. . 1000. 
*,800. ,600. ,500. ,4*0. , .00310, .00326, .00352, .00367, .00389, .00417, .00 
*435, .00456, .00495, .00538, .00585, .00663, .00791, .00898,4*0. , .00920, . 
*00955, .0101, .0106, .0112, .0123, .0133, .0147, .0173, .0202, .0231, .0274, 
*.0346, .0403,4*0. ,0.7021,20.94, 16*0.0) 

DATA ( ( (LR(N, J) , J-l, 18) ,N»1,4)= 

*6000. ,5000. ,4000. ,3000. ,2500. .2000. , 1500. , 1200. , 1000. ,800. .600. . 
*500. ,6*0. , .00308, .00310, .00309, .00309, .00322, .00352, .00420, .00491, 
* . 00562 , . 00662 , . 008 15,. 00930 , 6*0 . , . 00882 , . 00900 , . 00925 . . 00970 ..0104 
*0, .01205, .0151, .0182, .0215, .0264, .0343, .0405,6*0. .0.7273,20.65, 16* 
*0.0) 

DATA< < (MR(N, J) , J-l, 18) ,N-1,4>- 

*8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. , 1000. ,800. . 
*600 . , 500 . , 400 . , 4*0 . , . 00320 , . 00337 , . 00348 , . 00363 , . 00382 , . 00395 , . 004 
*10, .00443, .00497, .00567, .00672, .00834, .00960, .0113,4*0. , .00851, .00 
*900, .00931, .00972, .0104, .0112, .0123. .0142, .0167, . 0197, . 0242.  .  0314, 
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*.0372, .0457,4*0. ,0.7378, 18.57, 16*0.0) 

DATA  < ( <NR<N, J) , J-l, 18) ,N-1,4>- 

*10000. ,9000. ,8000. ,7000. ,6000. ,5000. ,4000. ,3000. ,2000. , 1500. . 1200. 
*, 1000. ,800. ,600. ,500. ,400. ,2*0. , .00295, .00299, .00303. .00310, .00318 
*, .00328, .00341, .00372, .00445, .00523. .00608, .00682. .00797, .00889. 
*.01101, .0129,2*0. , .00723, .00740, .00763, .00790, .00826, .00871. .00945 
*, .01085, .01370, .01645, .0195,-0228, .0278. .0357, .0419, .0511,2*0.  ,0.8 
*373,20.42, 16*0.0) 

DATA ( ( ( OR  <  N ,  J ) , J-l, 18) ,N=1,4)« 

*10000.  ,<9000.  ,8000.  ,7000.  ,6000.  ,5000.  ,4000.  ,3000.  ,2000.  ,  1500.  ,  1200. 
*, 1000. ,800. ,600. ,500. ,400. ,300. ,0. , .00294, .00302, .00309, .00317, .00 
*322,-00323, .00330, .00317, .00329, .00379, .00437, .00498, .00589, .00729 
*, .00833, .00980, .01215,0. . . 00716, . 00730, .00755, .00782, .00819, .00856 
*,  .00885,  .00956, .01145, . 01350, . 0159, .0181, .0220, .0285, .0336, .0411, 
*.0535,0. ,0.6541, 16.25, 16*0.0) 

DATA ( <  <PR<N, J) , j=l, 18) ,N=1,4)  = 

*8000. , 7000. , 6000. , 5000. , 4000. , 3000. , 2000. , 1500. , 1200. , 1000. , 800. , 
*600. ,500. ,400. ,300. ,200. .2*0. , .00302, .00312, .00322, .00333, .00344, . 
*00350, .00346, .00388, .00441, .00493, .00555, .00713, .00815. .00955. .011 
*85, .01600,2*0. , .00851, .00881, .00928, .00980, .01045, .01128, .01285, . 0 
*1475, .0170, .0195, .0238, .0306, .0359, .0437, .0566, .0811,2*0. .0.6541, 1 
*6.25, 16*0.0) 

DATA ( <  <QR<N,J) , J-l, 18)  ,N=1,4)  = 

*5000 . , 4000 . , 3000 . , 2000 . , 1 500 . , 1 200 . , 1 000 . , B00 . , 600 . , 500 . , 400 . , 300 . 
*,6*0. , .00281, .00281, .00263, .00268, .00294, .00338, .00379, .00448, .005 
*61, .00658, .00796, .01020,6*0. , .00809, .00835, .00875, .00962, .01088, . 0 
*125. .0144, .0178, .0232, .0275, .0339, .0442,6*0. ,0.6213, 13.27, 16*0.0) 
DATA ( ( < SR  <  N ,  J ) ,J  =  1, 18) , N= 1 , 4 ) = 

*3000. ,2000. , 1500. , 1200. , 1000. ,800. ,600. ,500. ,400. ,300. ,8*0. , .00277 
*, .00312, .00354, .00401, .00450, .00529, .00670, .00782, .00942, .01193, 
*8*0. , .00831, .00981, .01165, .0134. .0153, .0185. .0240, .0286, .0351,-046 
*0,8*0. ,0.6555, 13.80, 16*0.0) 

DATA< < <TR(N, J) , J-l, 18) ,N-1,4)- 

*3000. ,2000. , 1500. , 1200. , 1000. ,800. ,600. ,500. ,400. ,300. ,8*0. ,. 00293 
*, .00356, .00418, .00481, . 00545, . 00643, .00B02, .00922, .0110, .0138, 
*8*0. , .00981, .01185, .01395, .0162, .0189, .0230, .0302, .0361, .0448. .059 
*5,8*0. ,0.7352, 17.85, 16*0.0) 

DATA  <  <C<I)  ,  1-1,  14)=- 

*.375, .375, .500, .500, .375, .375, . 1875, .250. .250, .375, .375, .500, .750, 
*.750) 

DATA ( <D<I> , 1*1, 14)- 

*0.  ,  1.  ,0.  ,  1.  ,0.  ,  1.  ,0.  ,0.  ,2.  ,0.  ,2.  ,0.  ,0.  ,2.  ) 

DATA ( <S1 (I) , 1-1, 14)- 

*6.06,6.06,6.06,6.06,8.7,8.7, 11.1,11.1,11.1,11.1,11.1,11.1,11.1, 
*11.1) 

DATA  < ( <  VR  <  N , J ) ,J-1, 18) ,  N= 1 , 4 ) » 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. . 1000. 
*,800. ,600. ,500. ,4*0. , .00551, .00593, .00651, .00690. .00738, .00805, .00 
*849,  .00900,  .00970, .0104, .0112, .0124, .0144, .0160.4*0. , .0331, .0340, . 
*0354 , . 0363 , . 0375 , . 0394 , . 0406 , . 0426 , . 046 1,-0496,. 0532 , . 0587 , . 0682 . . 
*0755,4*0. , 1.2809,37.75, 16*0.0) 

DATA ( ( <WR<N, J) , J-l, 18) ,N-1,4)» 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. . 1500. , 1200. . 1000. 
*,800. ,600. ,5*0. , .00638, .00688, .00760, .00810, .00878, .00970, .0102. . 0 
*110, .0119, .0127, .0138, .0140, .0149,5*0. , .0494, .0510, .0531, .0547, .05 
*68, .0596, .0620, .0646, .0696, .0745, .0795, .0860, .0962,5*0. , 1.3119,57. 
*72, 16*0.0) 

DATA ( <  <XR<N,J) , J-l , 18) ,N-1,4)- 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. .2000. , 1500. , 1200.  ,  1000. 
* , 800 . , 600 . , 5*0 . , . 00568 , . 00605 , . 00655 , . 00690 , . 00734 . . 0079 1 , . 00829 , . 
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*00875, .00948, .0102, .0109, .0118, .0133,5*0. , .0300, . 0310, . 0322, . 0332, 
*.0347, . 0366, . 0381 , . 0402, . 0438, . 0474, . 0512, . 0571 , . 0667, 5*0. , 1 . 2700, 
*40.02, 16*0.0) 

DATA ( ( (YR(N, J) , J-l, 18) ,N»1,4)» 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. , 1000. 
*,800. ,600. ,5*0. , .00598, .00645, .00714, .00760. .00809, .00895, .00941, . 
*0100, .0108, .0113, .0118, .0122, .0128,5*0. , .0400, .0413, .0432, .0447. .0 
*463, .0491, .0511, .0540, .0588, .0634, .0680, .0752, .0880,5*0. , 1. 1270.56 
*.48, 16*0.0) 

DATA  < < (AS<N,J) , J= 1 , 18) ,N=1,4)  = 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. , 1000. 
*,800. ,600. ,500. ,4*0. , .00542, .00583, .00640, .00678, .00737. .00794, .00 
*835, .00885, .00951, .0103, .0112, .0126, .0149, .0169,4*0. , .0297,  .0306.  . 
*0319, .0328, .0340, .0359, .0374, .0394, .0430, .0472, .0515. .0585,  .0700,  . 
*0793,4*0. ,1.3108,40.83, 16*0.0) 

DATA  <  <  <BS(N,  J)  ,  J-l,  18)  ,IM=1,4)« 

*8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. . 1500. , 1200. , 1000. ,800. , 
*600 . , 500 . , 5*0 . , . 00630 . . 00690 , . 00730 ,. 00790 , . 00870 , . 00950 , . 00980 , . 0 
*106,. 0113,. 0121,. 0131, .0145, . 0154, 5*0. , .0340, .0395, .0410, .0428. .04 
*20, .0470, .0497, .0550, .0580, .0620, .0680, .0790, .0890,5*0. . 1.2351,45. 
*95, 16*0.0) 

DATA ( < (CS<N,J) ,J=1, 18) ,N=1,4)= 

*6000.  ,5000.  ,4000. .3000. ,2500. ,2000. , 1500. , 1200. , 1000. .800. ,600. . 
*500. ,6*0. , .00690, .00740, .00802, .00899, .00960, .0103, .0113, .0122, .01 
*30, .0142, .0161, .0177,6*0. , .0350. .0367, .0390, .0426, .0452, .0491, .055 
*3, .0610, .0662, .0738, .0848, .0925,6*0. , 1. 1946,44.25, 16*0.0) 

DATA ( ( (DS(N, J) , J-l, 18) ,N*1,4>= 

*8000. ,6000. ,5000. ,4000. ,3000. ,2500. .2000. , 1500. , 1200. , 1000. .800. , 
*600. ,500. ,5*0. , .00666, .00728, .00771, .00825, .00901, .00954, .0102, .01 
*12,-0119, .0125, .0137, .0155, .0168,5*0. , .0309, .0333, .0351, .0374, .040 
*8, .0461, .0464, .0512, .0558, .0600, .0670, .0772, .0850,5*0. , 1.3333,44.4 
*1, 16*0.0) 

DATA  <  (  (ES(N,J)  ,  J-l,  18)  ,IM-1.4)« 

*8000. , 6000. , 5000. , 4000. , 3000. , 2500. , 2000. , 1500. , 1200. , 1000. . 800. , 
*600. ,500. ,5*0. , .00701, .00761, .00800,-00853, .00922, .00972, .0103. .01 
*12, .0120, .0128, .0139, .0157, .0170,5*0. , .0349, .0364, .0375, .0390, .041 
*2, .0430, .0456, .0502, .0550, .0595. .0662, .0780, .0870,5*0. , 1.3332.44.3 
*9, 16*0.0) 

DATA ( < (FS(N.J) ,J«1, 18) ,N=1.4)= 

*10000. , 8000. , 6000. , 5000. , 4000. , 3000. , 2500. , 2000. , 1500. , 1200. , 1000. 
*,800. ,600. ,500. , 4*0. , . 00548, . 00588, .00645. .00684, .00793. .00811, .00 
*861, .00930, .0102, .0111, .0121, .0135, .0156. .0170,4*0. . .0242,  .0253, .0 
*271,  .0283,  .0300,  .0326. .0346, .0375, .0423, .0469, .0513, .0528. .0700, .0 
*796,4*0. , 1.3455,40.90, 16*0.0) 

DATA ( ( <GS(N, J) , J-l, 18) ,N-1,4)» 

*8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. , 1000. ,800. , 
*600. , 500. , 5*0. , . 00590. . 00650, . 00694, . 00752, . 00835, . 00889, . 00960.  . 0 
*105, .0112, .0119, .0130, .0148, .0161,5*0. , .0257, .0271, .0281, .0296, .03 
*19,  .0336,  .0363,  .0406,  .0442,  .0483,  .0550,  .0659,  .0741,5*0.  ,  l.?'TT4,3B. 
*30, 16*0.0) 

DATA  < <  <HS<N, J) , J-l, 18) ,N-1,4)- 

*8000. ,6000. ,6000. ,5000. ,4000. ,3000. ,2500. , 1500. , 1200. , 1000. .800. , 
*600. ,500. ,5*0. , .00557, .00604, .00640, .00680, .00739, .00777,  .00825,  .0 
*0B8B, .00950, .0104, .0117, . 0137, . 0150, 5*0. , .0220, .0233, .0242, .0255, . 
*0271,  .0283,  .0299,  .0332, .0368, .0410, .0474, .0570, .0641,5*0. , 1. 1843,3 
*3.125,16*0.0) 

DATA  < < (IS(N, J) , J-l, 18) ,N-1,4>- 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000.  , 1500.  ,  1200.  ,  1000. 
» , 800 . , 600 . . 500 . , 4*0 . , . 00432 , . 00462 , . 00508 . . 00537 , . 00576 . . 00630 , . 00 
*663, .00711, .00787, .00859, .00928, .0103, .0119, .0132,4*0. , .0151, .0158 
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*, .0170, .0178, .0190, .0208, .0222,-0244, .0289, .0314, .0370.  .0427,  .0516 
*, .0580,4*0. ,1.0299,29.00, 16*0.0) 

DATA ( ( (JS(N, J) , J  =  l, 18)  ,N»1,4)= 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. , 1000. 
*,800. ,600. ,500. ,4*0. , .00440, .00469, .00510, .00537, .00572, .00621, .00 
*655, .00699, .00762, .00831, .00894, .00981, .0112. .0122,4*0.  ,  .0156,  .016 
*8, .0175, .0183, .0194, .0213, .0227, .0248, .0288, .0313, .0362, .0416, .050 
*0,-0565,4*0. ,0.9455,28.25, 16*0.0) 

DATA ( <S2  < I ) , 1  =  1, 12)  = 

*11. 1, 12.22, 15.2, 13.95, 11.94, 15.4, 12. 18, 15.75,20.06, 19.82, 16. 12, 
*16.00) 

DATA ( <  <KS (N, J ) , J=l, 18) ,N=1,4)  = 

*8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. , 1000. ,800. , 
*600. ,500. ,5*0. , .00525, .00580, .00620, .00669, .00740, .00789,-00850, .0 
*0940,-0102, .0109, .0122, .0139, .0155,5*0. , .0197, .0209, .0218,  .0231,  .0 
*253, .0272, .0298, .0348, .0394, .0438, .0500, .0595,-0665,5*0. , 1.2012,33 
*.25, 16*0.0) 

DATA  (  (  < LS  <IM,J)  ,  J  =  1 ,  18)  ,  N=  1,4)  = 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. , 1000. 
*,800. ,600. ,500. ,4*0. , .00629, .00688, .00770, .00828, .00903, .0101, .010 
*8, .0119,-0133, .0146, .0156, .0171, .0192, .0205,4*0. , .0394, .0413, .0440 
*, .0458, .0487, .0530,-0560, .0607, .0680, .0752, .0826, .0942, . 113, . 130, 
*4*0. , 1.5887,65.00, 16*0.0) 

DATA  < <  <MS(N, J) , J=l, 18) ,N=1,4)  = 

*6000 . , 5000 . , 4000 . , 3000 . , 2500 . , 2000 . , 1 500 . , 1 200 , 1 000 . , 800 . , 600 . . 500 
*. ,400. ,300. ,4*0. , .00850, .00896, .00959. .01040, .01110, .01177. .01267, 
*.01327, .01373, .01427, .01520,-01580, .01675, .01810,4*0. . .0487, .0498, 
*.0516, .0540, .0558, .0584, .0628,-0676, .0726, .0800, .0913, . 1010, . 1145, 
*. 1390,4*0. ,0.9810,41.70, 16*0.0) 

DATA ( < (NS (N, J) , J=l, 18) ,N=1,4)= 

*6000.  ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. , 1000. ,800. ,600. . 
*500. ,400. ,5*0. , .01110, .01170, .01250, .0137, .0144, .0155, .0168, .0181, 
*.0192, .0204, .0223, .0233, .0247,5*0. ,. 0650, . 0664, . 0684, . 0712, . 0733, . 
*0765, .0817, .0870, .0927, . 1020, .1170, . 131, . 154,5*0. ,1.6373,61.60, 16* 
*0.0) 

DATA ( <  < OS  <N,J) , J= 1 , 18) ,N=1,4>  = 

*7000. ,6000. ,5000. ,4000. ,3000. ,2000. , 1500. , 1200. , 1000. , 800. , 600. , 
*500. , 400. , 300. , 4*0. , . 00452, . 00471 , . 00492, . 00522, . 00575, . 00682, . 007 
*44, .00830, .00911, .01045, .01255, .01415, .0166, .0205,4*0. , .0126,  .0131 
*, .0137,-0146, .0162,-0198, . 0231 , . 0265, . 0306, .0347,-0429, .0493,-0592 
*, .0758,4*0. , 1. 1110,22.74, 16*0.0) 

DATA ( < (PS(N, J) , J-l, 18) ,N=1,4)= 

*6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. , 1000. ,800. ,600. , 
*500. ,400. ,300. ,4*0. , . 00510, .00537, .00570, .00617, .00650, .00692,  .007 
*56, .00809,-00864, .00952, . 01 107, . 01227. .01407, .0169,4*0. , .0167, .017 
*5, .0186, .0202, .0213, .0228, .0255, .0283, .0314, .0362,-0443, .0507, .060 
*1,-0757,4*0. ,0.9159,22.74, 16*0.0) 

DATA ( <  <QS(N, j) , J=l, 18) ,N=1,4)  = 

*9000. ,8000. ,7000. ,6000. .5000. ,4000. ,3000. .2000. . 1500. , 1200. , 1000. . 
*800. ,600. ,500. ,400. ,300. ,2*0. , .00512, .00530, .00557, .00591, .00635, . 
*00692, .00782, .00933, . 01065, . 01 19, .0129, .0141, .0169, .0191, .0023. .02 
*78,2*0. , .0183, . 01B4, .0189, .0196, .0203, .0218. .0241, .0290. .0341, .038 
*8, .0438, .0490, .0592, .0695, .0808, . 1025,2*0. . 1.5067,29.60, 16*0.0) 
DATA ( < (RS(N, J) , J=l, 18) ,N»1,4>= 

*6000. ,5000. ,4000. ,3000. ,2000. , 1500. , 1200. , 1000. ,800. ,600. ,500. ,7*0 
*. , .00619, .00649, .00713, . 00813, . 00992, .01125, .0124, .0136, .0154, .018 
*5, .0209,7*0. , .0203, .0211, .0227, .0248, .0294, .0339, .0386, .0440, .0499 
*, .0608, .0690,7*0. , 1.6100,28.50, 16*0.0) 

DATA ( ( <TS(N,J) , J=1 , 18) , N=1 , 4) = 

*3000. ,2000. , 1500. , 1200. , 1000. ,800. ,600. .500. , 10*0. , . 00855. . 00995. . 
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*01115, .0120, .0129, .0144, .0173, .0197, 10*0. , .0309, .0349, .0387,  .0422, 
*.0459,  .0520, .0621, .0699, 10*0. , 1.5267,29.65, 16*0.0) 

DATA ( ( (US  <N,J) ,  J  =  1 , 18) , N=1 , 4) * 

*3000. ,2000. , 1500. , 1200. , 1000.  ,800. .600.  ,500.  .  10*0. ,  .00880,  .01015.  . 
*01155, .0128, .0139, .0154, .0180, .0202, 10*0. , .0420, .0450, .0492, .0535, 
*.0577, .0640, .0747, .0832, 10*0. , 1.5546,31.70, 16*0.0) 

DATA  <  (  <VS(N,  J)  ,  J=l.  18)  ,N=1,4)= 

*5000 . , 4000 . , 3000 . , 2500 . , 2000 .  , 1 500 . , 1 200 . , 1 000 . , 800 . , 600 . , 500 . , 
*400 . , 300 . , 5*0 . , . 0072 1 , . 00764 , . 00822 , . 00859 , . 00908 , . 00987 ,.01 060 , . 0 
*1123,  .01205, .01352, .0150, .0176, .0266,5*0. , .0310, .0315, .0334. .0357, 
*.0379, .0400, .0429, .0464, .0517, .0607, .0679, .0781, .0937.5*0. . 1.2249. 
*28. 11, 16*0.0) 

DATA  <  (  < WS  (N,  J)  ,J=1,  18)  ,  N=1 , 4)  =» 

*5000. ,4000. ,3000. ,2000. , 1500. , 1200. , 1000. ,800. ,600. ,500. ,8*0. . .007 
*78, .00838, .00925, .01085, .01205, .0132, .0142, .0159,  .0188,  .0209,8*0.  , 
*.0295, .0307, .0328, .0373, .0418, .0460, .0502, .0568, .0675, .0765,8*0. . 
*1 . 6197,33.45, 16*0.0) 

DATA ( <S3(I) ,  1  =  1,3)  = 

*11.44,11.5, 17.8) 

DATA< ( <YS(N, J) , J-l, 18) ,N=1,4)= 

*8000. ,6000. ,5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200.  ,  1000.  .800.  . 
*600. ,500. ,5*0. , .00712, .00794, .00846, .00920, .01025, .0110, .0119, .013 
*2, .0144, .0153, .0165, . 01 75, . 01 79, 5*0. . .0359, .0401, .0430. .0469, .0524 
*,  .0563,  .0615,  .0691, .0758, .0819, .0888, .0984, . 1045,5*0. , 1.3812,49.80 
*, 16*0. 0) 

DATA  <  (  <XS<IM,  J)  ,  J-l,  18)  ,N-1,4)» 

*10000. ,8000. ,6000. ,5000. ,4000. ,3000. ,2500. .2000. . 1500. , 1200. , 1000. 
*,800. ,600. ,500. ,400. ,3*0. , .00686, .00746, .00831, .00890, .00970, .0107 
*7, .01155, .0126, .0140, .0150, .0158, . 0167, . 0178, . 0185, .0194,3*0. . .033 
*1, .0357, .0398, .0427, .0467, . 0525. . 0567, . 0625, .0704, .0779, .0845, .092 
*6, . 1035, . Ill, . 118,3*0. ,1.2860,47.20, 16*0.0) 

DATA  (  <  (ZS(N,  J)  ,  J=l,  18),  N—  1,4)  = 

*5000. ,4000. ,3000. ,2500. ,2000. , 1500. , 1200. . 1000. ,800. ,600. ,8*0. , .00 
*675, .00740, .00835, .00900, .00982, .0110, .0120, .0129, .0142, .0158,8*0. 
*, .0293, .0320, .0358, .0385, .0421, .0478, .0530, .0579,  .0643,  .0738,8*0.  , 
*1.3911,44.28, 16*0.0) 

DATA ( < (UR(N, J) , J=l, 18) ,N-1,4>« 

*2000.  ,  1500.  ,  1200. , 1000. , 800. , 600. , 500. ,11*0.,. 00294, . 00349, . 00418, 
*.00482, .00581. .00735, .00856, 11*0. . .00118, .0135, .0157. .0183, .0228, . 
*0301, .0359, 11*0. ,0.6454, 17.85, 16*0.0) 

GD  TO  (100,1 000 , 3000 , 5000 ) , TYPE 

100  GO  TO  <110,120,130,140,150,200,210,220,230,240,250,260,270.280.290 
*,300,310,320) ,  NN 

110  CALL  INT <AR, RE, NST.F) .RETURNS (6000) 

120  CALL  INT (BR, RE, NST.F) .RETURNS (6000) 

130  CALL  INT (CR, RE, NST.F) .RETURNS (6000) 

140  CALL  INT (DR, RE, NST.F) .RETURNS (6000) 

150  CALL  INT (ER, RE, NST,F) .RETURNS (6000) 

200  CALL  INT <FR, RE, NST,F) .RETURNS (6000) 

210  CALL  INT ( IR, RE, NST, F) , RETURNS (6000) 

220  CALL  INT (JR, RE, NST, F) , RETURNS (6000) 

230  CALL  INT (KR, RE, NST.F) .RETURNS (6000) 

240  CALL  INT (LR, RE, NST, F) .RETURNS (6000) 

250  CALL  INT (MR, RE, NST, F) , RETURNS (6000) 

260  CALL  INT (NR, RE, NST, F) .RETURNS (6000) 

270  CALL  INT (OR, RE, NST.F) .RETURNS C6000) 

280  CALL  INT (PR, RE, NST.F) .RETURNS (6000) 

290  CALL  INT (QR, RE, NST, F) .RETURNS (6000) 

300  CALL  INT (SR, RE, NST.F) .RETURNS (6000) 

310  CALL  INT (TR, RE, NST.F) .RETURNS (6000) 


34 


PAS  82-41 


320  CALL  INT <UR, RE,  NST, F) .RETURNS (6000) 

1000  SO  TO  (1100,1200,1300,1400,1500,1600,1700,1800.1900,2000,2100, 
*2200, 2300, 2400) , NN 

1100  CALL  INT <VR, RE, NST, F) .RETURNS (6000) 

1200  CALL  I NT <WR, RE, NST,F) .RETURNS (6000) 

1300  CALL  INT ( XR, RE, NST, F) , RETURNS (6000) 

1400  CALL  INT (YR, RE, NST, F) .RETURNS (6000) 

1500  CALL  INT (AS, RE, NST, F) .RETURNS (6000) 

1600  CALL  INT (BS, RE, NST, F) .RETURNS (6000) 

1700  CALL  INT (CS, RE, NST, F) .RETURNS (6000) 

1900  CALL  INT (DS, RE, NST, F) .RETURNS (6000) 

1900  CALL  INT <ES, RE. NST, F) .RETURNS (6000) 

2000  CALL  INT (FS, RE, NST, F) .RETURNS (6000) 

2100  CALL  INT (SS, RE, NST, F) , RETURNS (6000) 

2200  CALL  INT (HS, RE, NST, F) .RETURNS (6000) 

2300  CALL  INT (IS, RE. NST, F) .RETURNS (6000) 

2400  CALL  INT <JS, RE, NST, F) , RETURNS (6000) 

3000  GO  TO  (3100,3200,3300,3400,3500,3600,3700,3800,3900,4000.4100,4200 
*)  ,  NN 

3100  CALL  INT  < KS , RE , N 3T , F ) , RETURNS  < 6000 ) 

3200  CALL  INT (LS, RE, NST, F) .RETURNS (6000) 

3300  CALL  INT (MS, RE. NST, F) .RETURNS (6000) 

3400  CALL  I NT ( NS , RE , NST , F ) , RETURNS ( 6000 ) 

3500  CALL  INT (OS, RE. NST, F) .RETURNS (6000) 

3600  CALL  INT (PS, RE. NST, F) .RETURNS (6000) 

3700  CALL  INT ( QS , RE . NST , F ) , RETURNS < 6000 ) 

3800  CALL  INT <RS, RE, NST, F) .RETURNS (6000) 

3900  CALL  INT <TS, RE, NST, F) .RETURNS (6000) 

4000  CALL  INT (US, RE. NST, F) .RETURNS (6000) 

4100  CALL  INT (VS, RE, NST, F) .RETURNS (6000) 

4200  CALL  INT (WS, RE, NST, F) .RETURNS (6000) 

5000  GO  TO  (5100,5200,5300,6000) ,NN 
5100  CALL  INT <YS, RE. NST, F) .RETURNS (6000) 

5200  CALL  INT <XS, RE. NST, F) .RETURNS (6000) 

5300  CALL  INT (ZS, RE, NST, F) .RETURNS (6000) 

6000  RETURN 
END 

*DECK  SURF 

SUBROUTINE  SURF  (TYPE, NS, A1 , B1 , SF, PS, RH, DEL, BET, FR, WF, WP> 

INTEGER  TYPE 
DIMENSION  S (4, 18, 12) 

DATA ( ( ( S ( 1 , I , N ) , N« 1, 12) , 1=1, 18)= 

*  1.,. 0000,0.,  2.00, .750, .04740, .032, 12.000,  76.1,-606,1.693,1.068 

*, 2. , .0000,0. ,  3.01, .750, .03546, .032, 12.000,  98.3, .706. 1.460, 1.107 

*, 3. , .0000,0. ,  3.97, .750, .02820, .032, 12.000,  119.4, .766, 1.353, 1. 145 
*, 4. , .0000,0. ,  5.30, .470, .02016, .006,  2.490,  188.0,-719,1.013,1.033 
*, 5. , .0000,0. ,  6.20, .405, .01820, .010,  1.200,  204. 0, . 728, 1 . 408, 1 . 066 
*, 6. , .0000,0. ,  9.03, .823, .01522, .008,  1.190,  244. 0, . 888, 1 . 010, 1 . 078 
*, 7. , .0000,0. ,11.10, .250, .01012, .006,  2.500,  367.0, .756, 1.032, 1.071 
*, 8. , .0000,0. ,11.11,. 480, .01 153, .008,  8.000,  312.0, .854, 1.000, 1.089 
* , 9. , .0000,0. , 14.77, .330, .00848, .006,  2.510,  420.0, .844, 1.018, 1.097 
1,0. , .0000,0. ,15.08,-418, .00876, .006,  6.840,  414. 0, . 870, 1 . 014, 1 . 100 

1. 1. .  .0000.0. .19.86. .250. .00615. .006,  2.510,  561.0, .849, 1.024, 1.135 

1.2. .  .0000.0. .10.27. .544. .01259. .010,  5.000,  289.9, .863, 1.018.1.114 

1.3. .  .0000.0. .11.94. .249. .00940. .006,  5.000,  393.0, .769, 1.023, 1.077 

1.4. .  .0000.0. . 12.00. .250. .00941. .006,  2.500,  392.7, .773, 1.023, 1.078 

1.5. .  .0000.0. . 16.96. .256. .00565,-006,  5.000,  607.8, .861, 1.009, 1.113 

1.6. .  .0000.0. .25.79. .204. .00377. .006,  2.500,  855.6, .884, 1.012. 1. 183 

1.7. .  .0000.0. .30.33. .345. .00401. .006,  2.500,  812.5, .928, 1.008, 1.223 

1.8. .  .0000.0. .46.45. . 100. .00264. .002,  2.630, 1332.5, .837, 1.020, 1. 102 
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*  ) 

DATA ( ( <S<2, I,N> ,N*1, 12) , 1-1, 14)« 

*  1.,. 3750,0.,  6.06, .250, .01460, .006,  0.055,  256. 0, . 640, 1 . 024, 1 . 038 
*,2. , .3750, 1. ,  6.06, .250, .01460, .006,  0.130,  256. 0, . 640, 1 . 024. 1 . 038 
*, 3. , .5000,0. ,  6.06, .250, .01460, .006,  0.055,  256. 0, . 640, 1 . 024, 1 . 038 
*,4. ,  .5000,  1.  ,  6.06, .250, .01460, .006,  0.130,  256.0,-640,1.024,1.038 
#, 5. , .3750,0. ,  8.70, .250,  .01196,  .006,  0.055,  307. 0,  .  705.,  1 . 024,  1 . 055 
*,6. , .3750, 1. ,  8.70, .250, .01196. .006,  0.080,  307.0,-705,1.024,1.055 
*, 7. , . 1875,0. ,11.10, .250, .01012, .006,  0.055,  367.0, .756, 1.024, 1.071 
*, 8. , .2500,0. ,11.10, .250, .01012. .006.  0.035,  367.0, .756, 1.024, 1.071 
*,9. , .2500,2. ,11. 10, .250, .01012, .006,  0.055,  367.0, .756, 1.024, 1.071 
1,0.,. 3750, 0. , 1 1 . 10, . 250, . 01012, . 006.  0. 055.  367.0, . 756, 1 . 024, 1.071 

1.1..  .3750.2. .11.10. .250. .01012. .006,  0.055,  367.0, .756, 1.024, 1.071 

1.2. .  .5000.0. .11.10. .250. .01012. .006,  0.055,  367.0, .756, 1.024, 1.071 

1.3. .  .7500.0. .11. 10. .250. .01012. .006,  0.040,  367.0, .756, 1.024, 1.071 

1.4. .  .7500.2. .11.10. .250. .01012. .006,  0.040,  367.0, .756, 1.024, 1.071 

*  ) 

DATA < <  <S<3, I,N> ,N=1, 12) , 1  =  1,  12)  = 

*  1. , .2500,0. , 11. 10, .250, .01012, .006,  0.250,  367. 0, . 756, 1 . 024, 1 . 07 1 
*, 2. , .0938,0. , 12.22, .485, .01120, .004,  0.094,  340.0, .862, 1.000, 1. 108 
*,3. , . 1250,0. , 15.20, .414, .00868,-006,  0.125,  417.0, .873, 1.015, 1.222 
*, 4. , . 1250,0. ,13.95,  .375,  .00879,-010,  0.125.  381.0, .840, 1.027, 1.385 
*,5. , .5000,0. ,11.94, .237, .00744, .006,  0.500,  461.0, .796, 1.512, 1.072 
*,6. , .2500,0. , 15.40, .206,-00527, .006,  0.250,  642.0, .816. 1.404, 1. 102 
*,7. , . 1667,0. , 12. 18, .353, .00885, .004,  0. 178,  422.0, .B47, 1.362, 1.051 
*, 8. , . 1429,0. , 15.75, .304, .00679, .004,  0. 143,  526.0, .859, 1.327, 1.067 
* , 9 . , . 1250,0. ,20.06, .201, .00489, .004,  0. 125,  69B.0, .843, 1.373, 1.087 
1,0. , . 1250,0. , 19.82, .205, .00509, .004,  0. 125,  680. 0, . 841 , 1 . 377, 1.086 

1.1.. . 1250.0. .16.12. .206. .00611. .006,  0.125,  660.0,-823, 1.381,1.107 

1.2. .  . 1250.0. . 16.00. .255. .00514. .006,  0. 125,  550.0, .845, 1.366, 1. 106 

*  ) 

DAT A<  <  <S<4, I , N ) ,  N=  1 , 12) , 1  =  1,  3)* 

*  1. , .3750,0. , 11.48, .413, .01060, .006,  0.078.  351 . 0, . 847, 1 . 015, 1 . 074 
*,2. , .3750,0. , 11.50, .375, .00993, .010,  0.078,  347.0, .822,1.027, 1. 130 
*, 3. , .3750,0. , 17.80, .413, .00696, .006,  0.078,  514.0, .892, 1.015,1.120 

*  ) 

A 1  =*S  <  TYPE ,  NS ,  2 ) 

B 1 *S  <  TYPE , NS , 3 ) 

SF — S <  TYPE , NS , 4 ) 

PS=S ( TYPE , NS , 5 ) 

RH-S  < TYPE , NS . 6 ) /4. 

DEL“S  <  TYPE , MS , 7 ) 

BET — S ( TYPE , NS , 9 ) 

FR=S (TYPE, NS, 10) 

WF*S (TYPE, NS, 1 1 ) 

WP*S ( TYPE , NS , 12) 

RETURN 

END 

*DECK  DEFAULT 

BLOCK  DATA  DEFAULT 

DEFAULT  VALUES  FOR  NAMELIST  INDATA 
TAFU  A 

COMMON / I NPUT / TYP A , NS A , TYPG , NSQ , RLEN I , RL I , NL . I AFRA , AFRA I , NA , W A , P I NA 
2  , TINA, PEXG,WB, TING, FAR, VINPUT 

DATA  TYP A, NSA, TYPG, NSQ  /l,  7,  1.  7  / 

DATA  RLEN I , RL I , NL , I AFRA , AFRA I , N A  /3.0,  1.0,  5,  25,  25.0,  4  / 

DATA  WA, PINA, TINA  /90.0,  116.4,  1040.5  / 

DATA  WG, PEXG, TING  /101.45,  14.90.  1646.4  / 
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DATA  FAR, VINPUT  /0.0145,  90.0  / 


END 

*DECK  INTERP 

SUBROUT I NE  I NT ( NR , RE , NST , F ) , RETURNS  <  A ) 

REAL  NR, NST 
DIMENSION  NR  < 4 , 18) 

N-0 

N*N+1 

J-0 

5  J»J+1 

I F ( RE . LT . NR ( N , J ) )  GO  TO  5 
IF ( J. EO. 1)  GO  TO  19 
I F ( NR  < 1 , J ) . EQ . 0 . )  GO  TO  10 
X*NR  < 1 , J  — 1 ) -NR ( 1 , J  > 

Y-NRI1,  J-l) -RE 
Z*Y/X 

NST«Z*<NR<2, J)-NR<2, J-l)  ) 

NST-NR<2,  J-D+NST 
F“Z*  <NR  <3, J) —NR (3, J-l ) ) 

F-NR<3, J-l)+F 
GO  TO  20 
10  F-NR(4,2)/RE 

NST-NR  <4, 1 ) /RE**0. 7 
GO  TO  20 

19  WRITE (6, 15)  RE 

15  FORMAT < IX, "REYNOLDS  NUMBER  OUT  OF  RANGE  OF  PROGRAMMED  TABLES  =  ", 
1F10. 1) 

NST-0. 

F»0. 

20  RETURN  A 
END 

*DECK  TRANSP 

SUBROUT I NE  TRANSP  <  T , FAR , CP , TK , MU , MW ) 

REAL  MU, MUA, MUF, MUG, MW 
C 

C******************#*******************************-********************* 

c 

DATA  A1 , A2, A3, A4, AS, A6, A7, A8/ 1 . 01 1 5540E-25 , - 1 . 4526770E-2 1 , 

A7 . 6215767E-18, -1 . 512B259E-14, -6. 7178376E-12. 6. 5519486E-08, 

B-5 . 1 536879E-05 , 2 . 502005 1 E-0 1 / 

DATA  B1 , B2, B3, B4, B5, B6, B7, B8/7. 267B710E-25, -1 . 333566BE-20, 

A1 . 0212913E— 16, -4. 2051 104E-13, 9. 9686793E-10, -1 . 3771901E-06, 

B1 . 2258630E-03, 7. 3B1663BE-02/ 

DATA  Cl , C2, C3, C4, C5, C6, 07/ -6. 2176401E-22, 7. 1827364E-18, 

A-3. 1410386E-14, 6. 7214720E-1 1 , -7. 5336781E-8, 6. 1979074E-5, 

B-4.81E-3/ 

DATA  D 1 , D2 , D3 , D4 , D5 , D6 , D7/1 . 0404582E-19, -7. 5213293E-16. 

A2. 1637607E-12, —3. 1593096E-9, 2. 4649233E-6, -9. 0B00204E-4, 1. 1073E-1/ 

DATA  El , E2, E3, E4, E5, E6, E7/2. 4724974E-21 , -1 . 6756272E- 17, 

A4. 1505396E-14, -3. 9906519E-1 1,-9. 1347177E-9, 8. B743855E-5, 

B2.98E-3/ 

DATA  F 1 , F2 , F3 , F4 , F5 , F6 , F7 / -2 . 0255 1 69E- 19,1.41 96996E- 1 5 , 

A-3. 9713025E-12, 5. 6582466E-9, -4. 3414613E-6, 1 . 8135009E-3, -3. 3929E-1/ 
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IF  <T  .LT.  500.  .OR.  T  .GT.  2000.)  GO  TO  100 
IF <FAR  -LT.  0.0  .OR.  FAR  .GT.  0.034826)  GO  TO  100 
CPA«<  <  (  <  <  <A1*T+A2)*T+A3)*T+A4)*T+A5)*T+A6)*T+A7)*T+A8 
CP=CPA 

IF (FAR. EQ. 0. )  GO  TO  30 

CPF= ( < ( ( ( <B1*T+B2) *T+B3) *T+B4) *T+B5) *T+B6) *T+B7) *T+B8 
CPG= <CPA+FAR*CPF) / ( 1 . +FAR) 

CP=CPG 

30  TKA=  <  < ( < ( C 1 *T+C2 ) *T+C3 ) *T+C4  >  *T+C5 ) *T+C6 ) *T+C7 
TK-TKA/3600. 

IF (FAR. EQ. 0. )  GO  TO  40 

TKF= ( ( ( ( <D1*T+D2) *T+D3) *T+D4> *T+D5> *T+D6> *T+D7 
TKG® <TKA+FAR*TKF) / ( 1 . +FAR) 

TK=TKG / 3600 . 

40  MUA» ( ( <  <  <E1*T+E2) *T+E3) *T+E4) #T+E5) *T+E6) *T+E7 
MU=MUA/3600. 

IF  <FAR. EQ. 0.0)  GO  TO  50 

MUF= ( ( ( ( <F1*T+F2) *T+F3> #T+F4) *T+F5>  *T+F6> *T+F7 
MUG=  <MUA+FAR*MUF> /  <  1 .  +FAR) 

MU=MUG / 3600 . 

50  MW“28. 97-. 946186*FAR 
RETURN 

100  WRITE  (6,  101)  T,  FAR 

101  FORMAT (10X.25HTRANSP  INPUT  OUT  OF  RANGE, 5X , 7HTEMP  =  , FB. 2, 5X , 6HFAR 
*  =  , F7. 4) 

RETURN 

END 

*DECK  BENDLOS 

SUBROUTINE  BENDLOS (X,Y) 

C 

Z=X*57. 29578 

Y- . 29227 1 3E-0 1 - . 2639695E-02*Z+ . 2272872E-03*Z**2- . 1 B50293E-05*Z  **3+ 
2  . 3655184E-07*Z**4-. 449784E-09*Z**5+. 208891  IE-1 1*Z**6 

C 

RETURN 

END 
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r; 

1  * 

- '  1 

APPENDIX  B 

j 

MAIN  PROGRAM  VARIABLES 

A 

plate  thickness,  In 

AA/G 

total  heat  transfer,  ft* 

ACA/G 

minimum  free  flow  area,  ft 

AFRA/G 

frontal  area,  ft2 

AFRAI 

value  of  frontal  area  Increase  per  Iteration,  ft2 

AHCMEXA/G 

header  core  matrix  exit  area,  ft^ 

AHCMINA/G 

header  core  matrix  Inlet  area,  ft2 

i 

ALHA/G 

ratio  of  total  transfer  on  one  side  to  total  volume  of  exchanger 

ANGEXA/G 

header  core  exit  angle,  radians 

ANGINA/G 

header  core  Inlet  angle,  radians 

AXA/G 

flns/ln 

BA/G 

plate  spacing.  In 

BETA/G 

BXA/G 

ratio  of  total  heat  transfer  area  to  volume  between  plates,  ft2/ft3 
denotes  different  fins 

i 

CA/G 

capacity  rate,  BTU/(hr  °R) 

t 

t 

CCA/G 

jet  contraction-area  ratio 

CFA/G 

pressure  effect  core  friction 

CMIN 

minimum  flow-stream  capacity  rate,  BTU/(hr  °R) 

COLBFA/G 

Colburn  factor 

CPA/G 

specific  heat,  BTU/{hr  °R) 

CR 

minimum  capacity  rate  ratio 

DELA/G 

fin  thickness.  In 

OELPA/G 

total  pressure  loss  on  one  side 

DELPAB/GB 

bend  pressure  loss 

DELPAB1/GB1 

Inlet  bend  pressure  loss  j 

DELPAB2/GB2 

exit  bend  pressure  loss 

DELPAC/GC 

core  pressure  loss 

DELPAH 

air-side  header  pressure  loss 

DELPT 

total  heat  exchanger  pressure  loss 

OEXA 

air-side  exit  header  diameter,  ft 

DINA 

air-side  Inlet  header  diameter,  ft  ! 

i 
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E 

ELA/G 

ERROR 

ETAFA 

ETAOA/G 

FA/G 

FAA/G 

FAR 

FLA 

FRA/G 

GA/G 

GC 

HA/G 

HFXA/G 

HINA 

I 

IAFRA 

ILA/G 

II 
K 

KA/G 

KBEXA/G 

KBINA/G 

KCA/G 

KDA/G 

KEA/G 

L 

LA/G 

LHGEXG 

LHGING 

LI 

MA/G 

MAL/GL 


heat  exchanger  effectiveness 

pressure  loss  exit  effects 

heat  exchanger  effectiveness  tolerance 

fin  effectiveness 

surface  effectiveness 

friction  factor 

pressure  effect  flow  acceleration 
fuel  to  air  ratio 

free  stream  counter-flow  length,  ft 
fin  area/total  area 

flow  stream  mass  velocity,  lbm/(hr  ft2) 

gravitational  constant  in  Newton's  second  law,  lbmft/(lb|r  sec2) 
unit  film  conductance,  BTU/(hr  ft2  °R) 

ratio  factor  used  in  core  pressure  loss  to  account  for  matrix 
headers 

air-side  dynamic  velocity,  lb^/ft2 
counter 

initial  frontal  area,  ft2 
pressure  loss  entrance  effects 
counter 

unit  thermal  conductivity,  BTU/(hr  ft2  ®R/ft) 
fluid  thermal  conductivity,  BTU/(hr  ft  °R) 
exit  bend  loss  coefficient 
Inlet  bend  loss  coefficient 
contraction-loss  coefficient 
momentum  velocity-distribution  coefficient 
expansion-loss  coefficient 
counter 

fin  length  from  root  to  center,  ft 
exit  gas-side  header  length,  ft 
Inlet  gas -side  header  length,  ft 
counter 

molecular  weight 

fin  effectiveness  parameter 
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MUA/G 

NA 

NARA/G 

NE 

NL 

NPRA/G 

NRA/G 

NSA/G 

NSTA/G 

NTU 

NXA/G 

OVALHT 

PAA/G 

PCDELPA/G 

PCDELPT 

PCNE 

PEXA/G 

PI 

PINA/G 

RA 

RHA/G 

RHO 

RHOEXA/G 

RHOINA/G 

RLEN 

RLENI 

RLI 

RU 

SFA/G 

SIGA/G 

SPVA/G 

SPVAM/GM 

TAVA/G 


viscosity  coefficient,  lbm/(hr  ft) 
number  of  frontal  area  Interatlons 
Reynolds  number  to  the  10"4 
heat  exchanger  effectlnvess 
number  of  length  Interatlon 
Prandtl  number 
Reynolds  number 

surface  number  used  as  Input  parameter 
Staton  number 

number  of  heat  transfer  units 
denotes  different  fins 

overall  length  of  heat  exchanger  enclosure,  ft 

percent  area  opening  on  one  side  to  total  area 

percent  total  pressure  drop  on  one  side,  % 

percent  total  pressure  drop  of  heat  exchanger,  % 

percent  heat  exchanger  effectiveness,  % 

outlet  pressure  on  one  side,  1 b^/1 n2 

*  -  constant  ■  3.1416 

Inlet  pressure  on  one  side,  lb^/in2 

parameter  used  in  overall  coefficient  of  heat  transfer, 

hr  ft2  °R/BTU 

hydraulic  radius,  ft 

density  of  heat  exchanger  material,  lbm/ft3 

3 

exit  density  of  fluid  on  one  side  of  heat  exchanger  headers,  lbm/ft 

Inlet  density  of  fluid  on  one  side  of  heat  exchanger  headers,  lbm/ft3 

flow  length,  ft 

Initial  flow  length,  ft 

value  of  length  Increase  per  Iteration,  ft 

gas  constant,  1n-lbf/(lbm  °R) 

flns/ln 

ratio  of  free-flow  area  to  frontal  area 
specific  volume,  ft3/lbm 
mean  conditions  specific  volume,  ft3/lbm 
average  temperature,  °R 
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TEXA/G  outlet  temperature,  °R 

TINA/G  inlet  temperature,  °R 

TLA/G  sum  of  pressure  loss  effects 

TYPA/G  fin  type 

UA  overall  coefficient  of  heat  transfer,  BTU/(hr  ft^  °R) 

VEXA  air-side  exit  header  velocity,  ft/sec 

VEXAC2/GC2  core  exit  velocity,  ft/sec 

VEXAH2/GH2  header  exit  velocity,  ft/sec 

VEXAM/GM  mean  exit  velocity,  ft/sec 

VINA  air-side  inlet  header  velocity,  ft/sec 

VINAC1/GC1  inlet  core  velocity,  ft/sec 

VINAHl/GHl  inlet  header  velocity,  ft/sec 

VINAM/GM  mean  inlet  velocity,  ft/sec 

VINPUT  input  velocity,  ft/sec 

VOL  volume,  ft^ 

WA/G  mass  flow  rate,  lb/sec 

WFA/G  factor  to  account  for  non-extended  fin  suface  weight 

WHXT  heat  exchanger  weight,  lb 

WIE  estimated  weight  of  header,  lb 

WPA/G  factor  to  account  for  plate  in  weight  calculations 

WPLA  total  weight  of  enclosure,  lb 

WTA/G  total  weight  of  fins  and  plate,  lb 

X  expontlal  variable  used  in  effectiveness  calculations 

XNCFL  heat  exchanger  non-flow  length,  ft 

NOTE: 

A/G  air  or  gas-side 
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APPENDIX  C 
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PAS  82-41 


DTNSRDC  PAS  82-41  ATTACHMENT  A 


KEY  TO  AIR-SIOE  OR  GAS-SIDE  FIN  GEOMETRIES* 

PLAIN  FINS  (TYPA/G  *  1) 

NSA/G  FINS/IN  PLATE  SPACING,  IN 


1 

2.00 

0.750 

2 

3.01 

0.750 

3 

3.97 

0.750 

4 

5.30 

0.470 

5 

6.20 

0.405 

6 

9.03 

0.823 

7 

11.10 

0.250 

8 

11.11 

0.480 

9 

14.77 

0.330 

10 

15.08 

0.418 

11 

19.86 

0.250 

12 

10.27 

0.544 

13 

11.94 

0.249 

14 

12.00 

0.250 

15 

16.96 

0.256 

16 

25.79 

0.204 

17 

30.33 

0.345 

18 

46.45 

0.100 

LOUVERED  FINS  (TYPA/G  =  2) 


NSA/G 

FINS/IN 

PLATE  SPACING,  IN 

LOUVER  SPACING, 

IN  LOUVER  GAP, 

1 

6.06 

0.250 

075 

0.055 

2 

6.06 

0.250 

0.375 

0.130 

3 

6.06 

0.250 

0.500 

0.055 

4 

6.06 

0.250 

0.500 

0.130 

5 

8.70 

0.250 

0.375 

0.055 

6 

8.70 

0.250 

0.375 

0.080 

7 

11.10 

0.250 

0.1875 

0.055 

8 

11.10 

0.250 

0.250 

0.035 

9 

11.10 

0.250 

0.250 

0.055 

10 

11.10 

0.250 

0.375 

0.055 

11 

11.10 

0.250 

0.375 

0.055 

12 

11.10 

0.250 

0.500 

0.055 

13 

11.10 

0.250 

0.750 

0.040 

14 

11.10 

0.250 

0.750 

0.040 

STRIP/OFFSET  FINS 

(TYPA/G  =  3) 

NSA/G 

FINS/IN 

PLATE 

SPACING, 

IN 

1 

11.10 

0.250 

2 

12.20 

0.485 

3 

15.20 

0.414 

4 

13.95 

0.375 

5 

11.94 

0.237 

6 

15.40 

0.206 

7 

12.18 

0.353 

8 

15.75 

0.304 

9 

20.06 

0.201 

10 

19.82 

0.205 

11 

16.12 

0.206 

12 

16.00 

0.255 

WAVY  FINS  (TYPA/G 

*  4) 

NSA/G 

FINS/IN 

PLATE 

SPACING, 

_IN 

1 

11 .44 

liRin 

2 

11.50 

0.375 

3 

17.80 

0.413 

*  From  Tables  9-3  (a),  (b),  (c)  and  (d)  In  Kay's  and  London  "Compact  Heat 
Exchangers." 


DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

V  DTNSRDC  REPORTS.  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT 

2.  DEPARTMENTAL  REPORTS.  A  SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY.  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION 

3.  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES.  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE  BY-CASE 
BASIS. 


